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Pre-face 

The  intent  of  this  study  was  to  take  a  -fresh  look  at  the  viability 
o-f  using  electric  thrusters  as  primary  propulsion  for  near-earth  space 
missions.  The  Shuttle  is  a  remarkable  system  and  enhancing  it  with  a 
reuseable  upper  stage  would  be  very  attractive. 

Chemical  rocket  systems,  even  with  SQA  LH  /  LOX  eny.nes  require 
enormous  amounts  of  fuel.  Without  aerobraking,  two  or  more  Shuttle 
refueling  flights  are  required  for  every  one  flight  bearing  GEO 
payload.  Upper  atmosphere  heat  transfer  and  drag  are  not  yet  well 
modeled,  so  the  technical  barriers  to  aerobraking  are  not  small. 
Electric  propulsion  makes  far  more  efficient  use  of  propellant,  but  has 
not  been  used  for  primary  propulsion  as  yet,  and  does  not  have  the 
production  base  that  chemical  has.  The  long  transfer  times  also  pose 
problems  for  revenue,  mission  promptness,  and  Van-Alien  degradation. 
Thus,  the  decisionmaker  has  no  easy  solution.  The  methodology 
developed  in  this  thesis  should  help  to  provide  information  as  to 
current  capabilities  of  optimized  ion  thruster  orbit  transfer  vehicles. 
It  also  should  permit  comparison  of  performance  and  cost  with  baseline 
chemical  rocket  systems  over  a  20  year  simulation  period. 

Originally,  I  had  hoped  to  include  self-field  magnetoplasmadynamic 
(MPD) ,  pulsed  inductive,  and  other  promising  thruster  technologies,  but 
will  have  to  leave  those  investigations  for  a  later  time  —  or  for 
future  classmates. 


I  heartily  acknowledge  the  assistance  and  advice,  both  technical 
and  paternal,  of  my  thesis  committee,  LTC  (Dr.)  Hark  Mekaru,  Dr. 


Wi  1 1  iam  Wiesel  ,  and  Or.  William  Elrod.  T^eir  enthusiasm  -for  the 
project,  encouragement,  and  humor  through  the  long  months  of  computer 
skirmishes  have  been  greatly  appreciated.  My  wife  and  2-yr  old  son 
have  been  extraordinarily  patient.  In  the  tradition  of  saving  the  best 
for  last,  I  wish  to  acknowledge  with  reverence  our  faithful  Savior,  the 
Lord  Jesus,  whose  wisdom  and  grace  are  matchless. 
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ABSTRACT 


'y  A  flexible  method&Togy  has  been ^ bevel  oped  for,-  optimizing  electric 
orbit  transfer  vehicles  (EOTVs)  and  comparing  them  with  baseline 
chemical  systems.  EOTVs  have  been  characterized  by  the  thruster 
technology  and  the  propellant  mass  versus  power  supply  mass  tor 
standardized  NASA  BIMQO  configurations.  Baseline  chemical  systems  are 
represented  by  the  Inertial  Upper  Stage.  X I US) , '  CENTAUR-G,  and  a 
proposed  reuseable  LOX-LH  Centaur  derivative. 

Five  electrostatic  propulsion  thrusters  were  chosen  for  the 
optimization.  These  were  the  baseline  NASA  /  Hughes  30-cm  J-Series 
Mercury  Ion  Thruster  and  four  derivatives.  Each  was  characterized 
through  linearization  of  experimental  data.  Relationships  of  input 
power  <KW)  to  the  thruster  vs  specific  impulse  and  input  power  vs 
thrust  were  developed.  The  first  relationship  along  with  equations  for 
power  supply  mass  and  propellant  mass  were  input  to  the  Sequential 
Unconstrained  Minimization  Technique  <SUMT)  nonlinear  ootimization 
program.  The  combination  of  the  propellant  mass  used  for  transfer  to 
GEO  and  return  and  the  power  supply  mass  was  minimized.  SUMT  runs  were 
made  for  the  five  thrusters  carrying  representative  payloads  from  1  to 
6  NavStar  GPS  satellites  with  associated  masses  of  90S  to  5443  KG. 
Transfer  times  were  then  calculated  for  each  of  these  payload  / 
optimized  ECTV  combinations.  Of  the  thrusters  chosen,  the  Ring-Cusp 
3-Grid  Xenon  thruster  accomplished  the  LEO  to  GEO  and  return  trios  with 
the  least  mass  and  the  minimum  transfer  time. 


Queueing-Graphical  Evaluation  Review  Technique  (QGERT)  program  was  used 


to  simulate  a  4-way  “f  1  y-off  “  between  EOTV,  I  US,  CENTAUR-G,  and  the 
Reuseable  Bi-Propellant  Vehicle  (RBPV)  .  The  results  of  the  29-year 
flyoff  comparison  were  used  to  assign  rough  Life-Cycle  Costs  (LCCs)  to 
the  operation  of  each  of  the  vehicles.  With  a  figure  of  merit  of  only 
LCC  of  each  system,  the  CENTAUR-G  appeared  best.  But,  when  using  a 
figure  of  merit  of  $LCC  per  KG  delivered  to  orbit,  the  EOTU  was  the 
best . 

Besides  the  initial  results,  the  methodology  can  be  used  by  those 
desiring  a  way  to  optimize  EGTVs  with  other  thruster  technologies  (eg. 
self-field  magn  e  topi  asmadyn  ami  c ,  pulsed  inductive).  Other  EQT'U 
configurations  may  be  used  for  the  optimizations  as  well  as  other 
payloads.  The  user  may  also  select  other  chemical  or  baseline 
systems  to  compare  with  EQTVs  in  the  QGERT  simulation. 
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ANALYSIS  OF  ELECTRIC  PROPULSION  ORBIT  TRANSFER  VEHICLES 


VS  I US,  CENTAUR-G, 

AND  A  REUSEABLE  B I  PROPELLANT  SYSTEM 


CHAPTER  I .  Background 


Introduction 

The  United  States  Space  Transportation  System  (STS)  is  an 
impressive  -feat  of  engineering  and  technology.  We  are  now  involved  in 
the  less  spectacular  business  of  practical  usage  and  operation  o f  the 
Shuttle.  It  has  been  made  clear  to  Headquarters  Air  Force  that 
enhancements  to  the  Shuttle  are  to  have  priority  when  considering 
future  space  systems.  In  keeping  with  this  charter,  this  thesis 
examines  the  improvement  of  Shuttle  capabilities  through  better 
transportation  in  the  near-earth  space  realm.  Expendable  upper  stages 
may  be  the  most  expeditious  and  the  only  realizable  way  to  do  business 
now,  but  will  continue  to  grow  unacceptably  expensive  as  more  numerous 
and  more  massive  space  systems  are  launched  and  placed  in  operational 
orbits.  Additionally,  the  capability  to  visit,  refurbish,  repair, 
investigate,  or  retrieve  space  assets  in  higher  orbits  is  non-existant 
as  of  this  writing. 

Very  briefly,  the  applicability  of  electric  propulsion  to  orbit 
transfer  missions  on  a  routine  basis  is  to  be  examined  during  mission 
simulations  of  20  years.  Optimal  parametric  trade-off  studies  will 
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first  oe  accomplished  to  define  the  eiectnc  propulsion  vehicle 
technologies  used  in  the  simulation.  The  performance  of  these 
electric  QTVs  (Orbit  Transfer  Vehicles)  will  then  be  compared  with 
current  baselines — the  operational  IUS  (Inertial  Upper  Stage),  the 
Centaur-G,  and  a  proposed  high-energy  liquid  bioropellant  OTV — in 
simulated  missions. 

The  results  of  the  trade-off  studies  and  the  simulations  will  be 
examined  in  light  of  enhanced  capability  for  the  Space  Transportation 
System  (STS)  and  for  reduction  in  the  number  of  Shuttle  launches  to 
place  given  satellite  constellations  and  space  structures  in  orbit. 


Hist  on 


Si tuation .  The  U.S.  Air  Force  has  been  "operational"  in  space 
for  a  number  of  years,  but  is  facing  a  new  era  with  the  advent  of  the 
Space  Shuttle.  A  recognized  weak  link  in  the  Space  Transportation 
System  (STS)  is  the  “upper  stage"  or  a  vehicle  to  place  satellites  in 
higher  orbits  than  the  Shuttle  orbiter  can  achieve.  Smaller 
satellites  sometimes  do  have  their  own  perigee  insertion  stage  and 

apogee  kick  motor  (AKM)  for  orbit  transfer.  But  satellites  are 
growing  larger  and  plans  are  being  made  for  Large  Space  Structures 
(LSS)  in  the  late  1980's,  1990's  and  beyond.  A  separate  Orbit 

Transfer  Vehicle  (OTV)  is  needed  not  only  for  these  larger  payloads, 
but  also  for  the  operational  capability  to  retrieve,  replace,  repair, 
or  refurbish  (R4)  satellites.  For  orbit  transfer,  the  Air 

Force-funded  Inertial  Upper  Stage  (IUS)  is  a  current  answer,  but  is 
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expensive  and  is  not  reuseable.  The  entire  vehicle,  including  costly 
3-axis  guidance  system,  is  discarded  after  one  mission  when  the  solid 
rocket  motors  are  spent. 

Much  more  efficient  would  be  a  vehicle  which  can  be  parked  in  LEO 
and  be  available  to  transport  Shuttle  payloads  to  (typically) 
Geosynchronous  Equitorial  Orbit. 


Importance.  The  need  for  a  reuseable  Orbit  Transfer  Vehicle  has 
been  highlighted  in  several  studies  (e.g.,  6,  9,  10,  16,  13,  19,  23, 

34,  35,  50,  51).  Current  solid-fuel  and  liquid-fuel  rocket 

technologies,  though  highly  advanced,  and  though  repeatedly  proven  in 
space,  require  a  significant  mass-fraction  of  propellant.  Solid  fuel 
motors  operate  typically  in  specific  impulse  ranges  of  approximately 
210-328  seconds  (24)  and  liquid-fuel  in  specific  impulse  ranges  of  390 
to  455  seconds  (numerous  Refs).  Hence,  the  mass  of  fuel  expended  to 
achieve  a  given  required  orbital  characteristic  velocity  increment, 
Av,  is  large.  The  "rocket  equation"  shown  below  illustrates  that  the 
relationship  between  exhaust  velocity  and  fuel  mass  is  exponential. 


Final  Mass 
_  -  Exp 

Initial  Mass 


Mission  Velocity 


Propellant  Velocity 


Initial  Mass  -  Final  Mass  s  Propellant  Mass 


(1) 


Electric  propulsion  is  attractive  in  space  because  it  makes  very 
efficient  use  of  fuel  which  allows  a  much  larger  payload  mass  to  be 
transported.  Also,  electric  propulsion  (EP)  typically  does  not 
generate  the  extremes  of  pressure  and  temperature  found  in  chemical 


propulsion.  (Temperature  is  an  exception  for  certain  types  of  EP) . 
Specific  impulses  may  range  from  1000  to  10,000+  seconds  +or  EP 


(numerous  Refs).  Electric  propulsion  also  has  advantages  of 
restar tabi 1 i ty  and  variation  of  thrust  level,  which  are  more  difficult 
with  chemical  propulsion.  Nuclear  propulsion  (i.e.  heating  of  H2  and 
expanding  through  a  nozzle)  still  has  some  of  the  limitations  of 
chemical  rockets  but  with  some  improvement  in  specific  impulse — to  300 
seconds. 

Thus,  electric  propulsion  has  the  edge  in  greatly  reduced  fuel 
mass  required.  And  added  to  this  inherent  advantage  is  that  of  having 
a  significant  amount  of  develoment — especially  for  ion -bombardment 
(e.g.,  12).  NASA-Lewis  Research  Center  has  developed  a  rather 
compl ete  modul ar  design  which  includes  the  engineering  detail  to 
submit  to  contractors.  But  the  problem  remains  that  an  EP  OTV  or 
vehicle  with  similar  capability  has  not  yet  been  produced.  Nor  has 
the  issue  of  an  appropriate  propulsion  system  been  resolved. 


Problem  Statement 

The  problem  is  to  pick  the  best  propulsion  technology  for  a 
reuseable,  modular  QTV.  Studies  to  date  have  concentrated  on 
trade-offs  to  see  which  technology  fits  each  mission  category  best. 
Typically  these  studies  have  generated  a  single  data  set  for  transfer 
time  or  a  graph  showing  technology  parameters  vs  v  or  payload  mass 
to  orbit.  The  issue  of  whether  electric  propulsion  can  do  the  orbit 
transfer  mission  has  been  well  studied  '.nearly  all  references),  but 
from  the  decision-maker's  standpoint,  the  question  of  which  EP 


needs 


more 


thruster  and  power  supply  -for  specific  missions 
clarification.  Several  optimum  parametric  OTV)  designs  CEP)  need  to  be 
determined  for  both  general  missions  and  for  specific  satellite 
constellation  emplacements  to  better  view  their  capability  from  the 
operational  standpoint. 

Thus,  the  problem  is  to  find  optimum  electric  propulsion 
technology  and  power  supply  technology  parameters  for  given  missions 
and  run  the  optimal  designs  in  18  to  40  year  simulations  to  examine 
long-haul  performance. 

Brief  Approach 

Orbital  parameters  and  payload  mass  for  projected  DOD  and  NASA 
missions  will  be  inputs  for  both  a  non-linear  optimization  program 
(3UMT) ,  and  a  simulation  program  CQGERT) .  The  SUMT  program  will 
determine  optimum  design  strategies  for  the  electric  OTVs  by  trading 
various  technical  and  performance  parameters.  The  simulation  runs 
will  in  essence  provide  a  parametric  "fly-off"  of  the  optimal  electric 
OTV)  and  the  IUS,  Centaur,  and  a  reuseable  bipropellant  system 
(projected).  Output  of  the  simulation  provides  number  of  satellites 
ferried  in  a  given  number  of  years,  average  wait  time  for  an  electric 
OTV)  in  LEO,  transfer  times,  and  optimum  number  of  vehicles  to  do  the 
mission.  This  will  form  a  basis  for  cost  comparison  between  OTVs  and 
stages.  The  results  should  give  the  decision-maker  a  good  handle  on 
technologies  to  invest  in  and  how  such  systems  should  perform  over  the 
1 ong  haul . 
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Scope  of  the  Literature  Search  and  Data  Base 

The  literature  search  has  been  quite  extensive  and  it  will, 
hopefully,  prove  to  be  a  significant  and  valuable  reference  for  the 
serious  reader.  However,  it  cannot  presume  to  be  exhaustive  as  there 
are,  undoubtedly,  many  fine  reports  and  studies  which  have  not  come  to 
the  attention  of  this  author.  Another  factor  which  affects 
comprehensiveness  is  the  on-going  electric  propulsion  research  at 
various  laboratories,  universities,  and  aerospace  firms.  Besides 
basic  research,  this  author  is  aware  of  continuing  studies  and  mission 
analyses  nearing  completion,  eg.,  Ref.  34.  Though  necessitating 
caveats,  such  continuing  work  tends  to  show  the  validity  of  this 
research  area.  The  sheer  number  of  contract  reports  and 
academic/professional  papers  show  electric  propulsion  to  be  a  viable 
area  of  space  propulsion.  It  is  difficult  to  ignore  such  efficient 
use  of  propellant  mass. 

The  following  list  delineates  areas  specifically  not  treated 
in  the  thesis  and  the  literature  search: 

1.  Detailed  derivation  of  the  theoretical  basis  for  electric 
propulsion  techniques. 

2.  Detailed  design  aspects  of  cryogenic/solid/nuclear/photon 
propulsion  systems. 


3.  Interpl anetary/interstel 1 ar  missions. 

4.  Every  concept  and  technique  of  electric  propulsion,  onlv 


those  applicable  to  the  given  mission  are  treated. 

5.  Detailed  analytical  derivation  of  orbit  transfers. 

6.  Data  on  materials  science  experiments  with  electric 
propulsion  components. 

7.  Complex,  stochastic  cost  estimation  techniques. 

3.  Detailed  reliability  forecasting. 

9.  Intricate  design  trade-offs  for  power  supplies — solar, 
nuclear,  thermionic,  fuel  cell,  battery,  etc. 

Solar  shadowing,  Beginning  Of  Life  (BOL)  and  End  Of  Life  (EGL) 
tradeoffs  for  solar  panels  are  also  left  to  other  studies,  al thougn 
sensitivity  to  these  effects  has  been  investigated  in  Appendix  V. 

Additionally,  the  OTV  “designs"  are  not  rigorous  engineering 
designs  ready  for  the  draftsman,  but  are  conceptual.  They  exist  in 
terms  of  the  parameters  of  the  system  (such  as  overall  weight, 
specific  impulse,  propellant  characteristics,  power  plant  efficiency, 
propulsion  efficiency,  thrust,  and  lifetime). 


MAibQd  of  Treatment  and  Organization 

All  58  of  the  bibliographical  entries  for  this  review  have  been 
grouped  under  six  headings.  Mi  thin  these  headings  are  sub-headings 
which  roughly  follow  the  development  of  the  thesis.  Not  every  one  of 
the  53  entries  will  be  reviewed,  since  several  are  general  references, 
such  as  textbooks,  since  others  are  ancillary  with  regard  to  the  main 
thrust  of  the  thesis,  and  since  several  references,  though  very 
relevant,  were  received  just  prior  to  publication.  Also,  conciseness 


is  a  courtesy  to  the  reader  as  well  as  a  stated  requirement. 

Therefore,  to  help  guide  you  in  the  literature  discussion,  the 
following  outline  is  given  <references  in  brackets  are  the  ones  not 
discussed} : 

1.  Background  Sources  <22,  23,  47,  48,  56} . 

History  (24) 

Advantages  of  Electric  Propulsion  (24) 

2.  Generic  Electric  Propulsion  <EP)  Sources  <3,52,14, 
17,31,32,38,33, 48} . 

Electric  Propulsion — A  Mature  Technology  (51) 

3.  Studies  of  Electric  Propulsion  as  Applied  to  Orbit 
Transfer  <7,3, 15, 27, 23, 43, 44, 54>. 

Chronological  Presentation  (13,9,6,35,18,19,50,34) 

4.  Data  Base  Sources  <5, 11, 13, 21, 30, 39, 40, 41, 42, 45,49, 55>. 

Overall  Design  Approach  Data  (10) 

Input  for  the  Trade-Off  Studies  (12) 

Input  for  the  Simulation  (29) 

5.  Sources  for  Analytical  Techniques  —  Multi-Criteria 
Trade-Offs,  Non-Linear  Optimizations,  (S'JMT)  ,  and 
Simulation  ( QGERT)  <4,58,37,20 ,33, 57> . 

8.  Orbital  Dynamics  Considerations  <2,25,53}. 

Low  Thrust  Transfer  is  Workable  (1) 

Please  note  that  the  references  are  not  locked  in  to  these  headings. 
For  instance,  reference  51  also  contains  excellent  background 
material.  Thus,  these  headings  serve  only  as  a  convenient  means  for 
organizing  the  entries  and  are  not  titles  for  thesis  chapters. 


Background  sources 

There  must  be  a  genuine  need  or  an  established  requirement  to 
justify  the  proposal  of  a  new  vehicle  or  system.  To  lay  a  foundation 
for  the  thesis,  the  fundamental  advantages  of  electric  propulsion  for 
space  missions  must  first  be  made  clear. 

History.  Dr.  Harold  Kaufman  describes  the  "Origin  of  the 
Electron -Bombardment  Ion  Thruster"  (26)  during  the  late  1958's.  This 
particular  thruster  is  relevant  to  this  thesis  in  that  it  has  become 
the  most  developed  form  of  electric  propulsion.  This  type  thruster 
now  bears  his  name  in  most  current  reports.  His  article  is  a 

historical  development  of  the  very  first  operable  el ectr on -bombardment 
ion  thruster  which  he  invented  or  developed  along  with  Milliam 
Kerslake  and  other  colleagues  while  at  NASA-Lewis  Research  Center  in 
the  late  1950s  and  early  1960s.  The  thruster  was  first  test-fired  in 
1960.  The  article  gives  an  appreciation  for  the  ingenuity  exercised 
in  its  crude  beginnings,  and  gives  an  appreciation  for  the  degree  of 
maturation  since  then  in  this  technology.  The  30-cm  electron 
bombardment  thruster  (12)  and  its  derivatives  are  considered 

state-of-the-art  (SQA)  as  of  this  writing. 


Advantages  of  Electric  Propulsion.  Dr.  Robert  6.  Jahn,  in  his 
textbook  used  at  Princeton,  "Physics  of  Electric  Propulsion,"  (24) 
provides  a  more  extensive  discussion  of  the  advantages  of  electric 
propulsion.  The  text  provides  the  theoretical  physics  background  for 
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electric  propulsion  and  the  electrical  acceleration  of  gases.  It  also 
provides  a  good  review  of  the  types  of  thrusters  and  their  techniques 
of  acceleration.  Chapter  1  is  referenced  as  background  for  this 
thesis  because  it  specifically  shows  the  "province  of  electric 
propulsion"  to  be  "high  impulse  space  missions."  It  explains  that  the 
primary  advantage  enjoyed  by  EP  is  a  much  higher  specific  impulse  than 
is  possible  with  standard  chemical  propulsion.  This  directly  impacts 
the  amount  of  propellant  necessary  to  carry  out  a  given  space  mission. 
Another  very  important  part  of  this  chapter  for  this  thesis  is  its 
discussion  of  the  power  supply  penalty.  That  is,  the  tradeoff  of 
power  supply  weight  for  higher  thrust  and  also  for  higher  specific 
impulse.  It  introduces  the  parameters  of  specific  power-plant  mass, 
conversion  efficiency,  and  character istic  velocity  increment,  Av. 
Also  discussed  are  some  specific  types  of  missions  where  electric 
propulsion  is  the  logical  choice  over  chemical  rockets. 

Generic  Electric  Propulsion  Sources 

The  majority  of  the  bibliographical  entries  for  this  thesis 
contains  information  regarding  types  and  classes  of  EP.  Of  interest 
to  this  thesis  are  those  classes  and  specific  thruster  designs  which 
are  most  promising  in  the  near  term  for  accomplishing  space  missions. 

Electric  Propulsion — A  Mature  Technoloov.  In  his  article, 
"Electric  Propulsion  Ready  for  Space  Missions",  (51),  Or.  Stuhlinger 
gives  an  excellent  overview  of  the  mature  state  of  many  EP  systems  and 
includes  many  details  of  performance.  Areas  discussed  included, 
"overview  of  EP  programs*  and  details  of  work  in  Japan,  West  Germany, 
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Great  Britain,  and  the  USSR,  as  well  as  U.S.  programs.  Organizations 
and  companies  with  EP  programs  included:  US  Air  Force,  NASA, 
Fairchild,  Phrasor  Technology,  Technion,  Hughes,  Lockheed,  and  several 
universities.  Other  areas  discussed  were,  "Spacecraft  Systems,” 
"Propulsion  Concepts,"  “Thruster  Design  and  Analysis,"  "Thruster 
Performance  and  Qualification,"  "Endurance  Tests,"  "Component  R&D," 
etc.,  all  of  which  point  directly  to  the  maturity  of  the  technology. 

Perhaps  the  most  applicable  part  of  Stuhlinger's  article  is  the 
extensive  treatment  and  discussion  of  the  30-cm  Kaufman  thruster 
because  it  is  the  most  developed  and  is  closest  to  being  ready  for 
orbit  transfer  missions. 


rich  in  material  and  data  for  the  thesis,  and  because  each  is  a  rather 
complete  analysis  in  itself,  the  following  discussion  will  primarily 
address  their  conclusions.  For  the  reader's  benefit,  attempt  has  been 
made  to  stick  to  the  essentials. 

D.G.  Fearn  (18)  found  that  two  uses  of  EP  were  especially 
attractive.  For  satellites  under  I860  kg,  the  EP  system  could  be  an 
integral  part  of  the  spacecraft  owing  to  the  small  fractional  mass  of 
fuel  required.  The  other  use  was  a  separate  OTV  or  tug  with  solar 
arrays  for  power.  He  pointed  out  the  possible  problem  of  solar  panel 
degradation  in  the  Uan  Allen  Belts,  but  still  found  the  application 


O.C.  Byers.  et  ai  ,  »9>  presented  *  general  -netnoaol  ogv  tor 
predicting  the  overall  EP  system  properties,  sucn  as  input  Dower  ar.o 
mass,  when  mission  parameters  and  propellant  type  were  given.  This 
provides  information  ■for  the  trade  studies  portion  of  this  thesis.  It 
aids  in  the  selection  of  the  right  vehicle  parameters. 

D.C.  Byers,  in  a  subsequent  study  16) ,  specifically  addressed  the 
major  theme  of  this  thesis,  "Upper  Stages  Utilizing  Electric 
Propulsion."  He  used  the  methodology  established  previously  (9)  to 
define  the  electric  thrust  system  and  its  power  requirements  in 
detail.  With  that  he  presented  the  payload  capabilities  of  upper 
stages  using  EP  for  LEO  to  GEO  orbit  transfer  missions.  This  thesis 
draws  upon  the  payload  data  from  this  paper  and,  additionally, 
accomplishes  the  simulation  of  an  operational  OTV  system  for  a  20  year 
period,  which  was  not  part  of  Byers'  paper. 

William  Pipes  headed  the  Martin  Marietta  Corp.  team  which  finished 
an  extensive  contract  study  (35)  for  the  AF  Rocket  Propulsion 
Laboratory  in  July,  1930.  Unlike  the  3  studies  discussed  so  far,  tnis 
one  considered  liquid  propellant  (both  storable  and  cryogenic),  and 
solid  propellant  as  well  as  electric  propulsion  systems  for  orbit 
transfer.  It  presented  the  results  in  terms  of  the  relative 
advantages  of  each  system  for  certain  weight  classes  of  satellites. 
Highl ighted  was  the  economic  benefit  of  Magnetopl asmadynamic  (MPD) 
propulsion  over  ion-bombardment,  although  significant  technical 
difficulties  (e.g. , electrode  erosion)  remain  with  MPD  systems. 

The  second  technical  report  by  D.G.  Fearn  (16)  tends  to  parallel 
the  Martin  Marietta  study  in  its  extensive  treatment  of  types  of 
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electric  propulsion  to  be  considered  for  orbit  transfer  missions.  The 
clear  conclusion  is  that  the  high  specific  impulse  offered  by  EP 
provides  an  enormous  economic  advantage  over  chemical  rocket  systems 
for  the  movement  of  non-priority  cargo.  Non-priority  is  specified 
because  the  transfer  times  are  typically  100  days  or  more.  For- 
projecting  performance  for  larger  diameter  ion  thrusters,  the  tables 
on  page  1 6  (16)  should  be  helpful. 

Robert  Finke  (19)  has  edited  a  volume  of  collected  papers  on  EP 
from  the  1979  AIAA  International  Electric  Propulsion  Conference.  This 
represents  one  of  the  most  complete  and  current  collections  available 
on  EP,  as  of  this  writing.  Its  comprehensiveness  is  probably  the  most 
valuable  contribution  to  this  thesis.  Several  papers  dealt  with 
application  of  EP  to  orbit  transfer  missions  and  especially  useful 
were  two  papers  dealing  with  cost-effectiveness. 

Capt  Jess  Sponable  (50)  has  shown  through  computer  modeling 
techniques  that  the  Space  Transpor tat  ion  System  (STS)  can  be  enhanced 
and  optimized  for  NASA/DOD  missions  through  the  employment  of  a  LEO 
space  station  and  OTV.  Although  not  the  primary  issue  in  his  thesis, 
he  also  demonstrated  significant  reductions  in  launch  rates  if  a  space 
station  and  an  EP  OTV  are  both  used.  My  thesis  will  not  include  a 
space  station  as  part  of  the  scenario  and  will  attempt  to  answer 
whether  similar  reductions  in  Shuttle  launch  rates  are  possible  with 
the  EP  OTV  alone.  Also  useful  from  his  thesis  was  the  NA3A/DGD 
mission  model - a  starting  point  for  developing  my  simulation  input. 

Capt  David  Perkins  is  in  the  process  of  having  his  report  printed 
which  is  entitled,  "Preliminary  Analysis  and  Comparison  of  Recoverable 


Space  Baaed  Orbit  Transfer  Vehicles  -tor  LEO  to  GEO  Missions."  (34) 
Publication  ia  forthcoming  in  late  1983  and  ahould  be  available  to 
certain  00D  users  at  that  time.  Needless  to  say,  this  report  ia  one 
o f  the  moat  current.  It  offers  a  look  at  aome  very  new  technologiea 
in  the  realma  of  power  suppl iea  and  thruster  concepts  rivaling  EP  in 
eff iciency .  While  this  report  does  not  perform  a  simulation  as  in 
thia  thesia,  it  provides  some  excellent  data  with  which  to  compare 
trade  studies. 


Data  Base  Sources 

Sources  to  this  point  have  certainly  contained  pertinent  data,  but 
the  following  sources  have  especially  useful  data  as  outlined  below. 

Overal 1  Design  Approach  Data.  Cake,  et  al ,  (18)  presented  a 
modular  approach  to  designing  an  EP  space  vehicle  which  used  solar 
arrays  for  power.  The  study  was  not  limited  to  an  orbit  transfer 
vehicle  but  was  flexible  enough  to  include  interplanetary  missions. 
Data  useful  for  this  thesis  includes  the  component  and  subsystem 
arrangements  necessary  for  an  operational  vehicle.  The  report 

emphasized  the  structural  and  thermal  integration  of  modular 

subsystems.  Three  approaches  to  a  Solar  Electric  Propulsion  (SEP) 
module  were  compared  on  the  basis  of  mass,  cost,  testing,  interfaces 
with  spacecraft,  simplicity,  maintainability,  and  reliability.  All 
portions  of  the  generated  data  have  relevence  to  my  data  base. 

Input  for  Trade  Studies.  The  "38-Centimeter  Ion  Thrust 
Subsystem  Design  Manual”  (12)  contains  specific  details  of  thrust 
output,  wattage  input,  specific  impulses  achieved,  and  total 


-  14  _ 


.*•  .N 


V 


L 
23 


eff iciencv.  Because  this  manual  was  prepared  with  the  goal  of 
producing  a  worKing  vehicle  with  necessary  specifications  tor 
contract,  it  has  data  for  several  levels  of  detail.  Drawings,  mass 
data,  configuration  layouts,  and  engineering  details  are  used  directlv 
in  this  thesis  to  characterize  the  Electric  Orbit  Transfer  vehicle 
(EOTV) .  As  stated  earlier,  this  is  the  most  highly  developed  of  the 
EP  thrusters  and,  hence,  one  of  the  most  likely  candidates  for  the 
first  generation  EOTV. 

Input  for  the  Simulation.  Ker slake  (2?)  has  covered 

performance  and  durability  of  an  ion  thruster  system  which  has  been  in 
orbit  and  operating  for  11  years  until  it  ran  out  of  fuel  in  1931. 
This  experiment  was  called  3ERT  II  (Space  Electric  Rocket  Test  -  II). 
This  report  serves  as  a  major  indicator  for  thruster  /  EP  system 
lifetimes  in  the  simulation  portion  of  the  thesis.  Eleven  years  is  a 
long  time  to  endure  the  extremes  of  space  and  this  test  provides 
empirical  evidence  and  credibility  for  the  concept  of  an  electric 
propulsion  GTV  in  addition  to  valuable  data  for  modeling. 


Orbital  Dynamics  Considerations 

Low  Thrust  Orbit  is  Workable.  Capt  Salvatore  Alfano  (1)  has 


shown  in  his  recent 

analytical  thesis  that 

the 

1  ow 

thrust 

orbit 

transfer  spiral  can 

be  solved  in  closed  form. 

He 

has 

verified 

the 

resulting  equations 

using  numerical  computer 

techniques 

.  Also, 

the 

results  were  compared  with  the  standard  Hohmann  transfer  ellipse.  His 
conclusion  was  that  the  resulting  spiral  transfer  is  optimal  for  low 
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thrust,  continuous  thrusting  and  thus  accomplishes  the  mission  in 
minimum  time  using  minimum  -fuel.  This  is  the  type  of  orbit  transfer 
that  must  be  used  with  electric  propulsion. 

Literature  Review  Conclusion 

The  literature  review  began  with  the  early  Kaufman  tnruster  and 
discussed  the  advantages  of  higher  specific  impulse  available  with 
Electric  Propulsion  (EP)  over  chemical  rocKets.  Then,  Stuhl  inger's 
article  was  reviewed,  pointing  out  that  EP  is  ready  for  implementation 
in  space  vehicles.  Many  studies  applying  EP  to  the  orbit  transfer 
mission  were  then  discussed  in  chronological  order  (1978  to  1983). 
Data  base  sources  were  discussed  last.  For  each  report  reviewed,  some 
indication  of  its  content  and  scope  was  given  as  well  as  its 
applicability  to  this  present  thesis. 

In  conclusion,  the  reader  should  realize  by  now  that  electric 
propulsion  has  received  a  significant  amount  of  attention  in  the 
scientific  community.  The  reasons  center  around  its  potential  for 
space  missions  requiring  high  specific  impulse  and/or  long  duration 
use.  Finally,  the  information  and  data  represented  in  this 
bibliography  provide  a  firm  foundation  for  the  thesis. 
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CHAPTER  III.  Approach  /  Methodoioc 


Research  Questions 

The  two  most  general  questions  to  be  answered  are,  "Which  electric 
propulsion  technology  is  the  best  choice  -tor  a  reuseable  OPJ,"  and 
"How  does  this  electric  system  compare  with  baseline  chemical  systems 
in  mission  performance  and  cost?"  The  more  specific  questions  to  be 
answered  are  as  follows: 

1.  Which  electric  thruster  technology  among  several  lab 
demonstration  prototypes  would  optimize  an  OTV  in  terms  of  reduced 
fuel  mass  and  reduced  power  supply  mass  for  given  missions? 

2.  Given  specific  missions,  what  are  the  transfer  times  and 
round-trip  mission  times  for  the  optimized  electric  OTVs? 

3.  Does  one  electric  thruster  technology  clearly  outperform  al 1 
others  for  each  mission? 

4.  Using  a  comparison  "fly-off"  simulation  between  Electric  OTVs 
(EOTVs)  and  baseline  chemical  systems  (IUS,  Centaur,  Reuseable 
8i -propel  1  ant) ,  is  a  reduction  in  shuttle  launch  rate  possible  using 
EOTVs? 

5.  Using  the  same  fly-off,  is  any  reduction  in  cost  over  present 
operations  suggested? 

6.  What  would  an  optimal  number  of  EOTVs  and  Re-useable  chemical 
OTVs  be  —  i.e.,  fleet  size? 

7.  Are  Shuttle  enhancements,  expanded  near-earth  operations,  or 
new  DOD  missions  suggested  by  the  performance  of  optimizeo  EOTVs? 


-  17  - 


wvv"  #7»r  v 


In  answering  these  questions,  the  algorithms  and  methodology  which 
are  to  be  developed  should  prove  helpful  to  decision-makers  and  to 
staff/commi  t  tees  considering  the  best  way  to  exploit  near-earth  space 
on  a  limited  budget.  A  set  of  research  objectives  has  been 
established  to  further  define  and  specify  how  these  questions  are  to 
be  answered. 


Research  Objectives 

The  following  research  objectives  clarify  what  this  tnesis  effort 
is  to  accomplish  in  answering  the  questions  posed  in  the  previous 
section.  This  should  also  help  to  narrow  the  scope  and  further  define 
the  problem. 

1.  Choose  at  least  4  electric  thruster  technologies  that  have 
been  operated  for  lab  testing  and  represent  developable  systems  in  the 
near  term  (5-15  years),  then  optimize  input  power  <kw>  and  specific 
impulse  (sec)  to  minimize  both  propellant  mass  and  power  supply  mass 
for  EOT'  3  using  these  thrusters  to  accomplish  given  missions. 

2.  Develop  an  algorithm  /  program  to  find  optimal  specific 
impulses  and  then  use  the  results  to  calculate  the  transfer  times  out 
to  an  orbit,  back  to  the  original  or  another  orbit,  and  the  combined 
round-trip  time.  This  should  be  flexible  enough  to  use  for  other 
missions  besides  deployment  of  satellites. 

3.  Examine  the  results  of  the  optimization  in  trading-off 
different  thrusters,  overall  EOTV  masses,  and  transfer  times  for  all 
given  missions  to  determine  if  one  technology  clearly  is  best. 


4.  Develop  a  GGERT  simulation  program  for  each  type  of  Upper 


Stag*  /'  OTV  and  use  it  to  fina  the  number  of  Shuttle  launches  required 
to  place  given  missions  in  their  -final  orbits. 

5.  Using  the  simulation  program,  determine  an  optimum  -fleet  size 
for  both  the  EOTV  and  the  Reuseable  Bi-Propel 1  ant  Vehicle  <RBPV) . 
“Optimum*  is  that  minimum  number  which  will  preclude  satellites  / 
payloads  from  queuing  in  low  earth  orbit  while  waiting  for  transfer, 
and  provide  some  vehicle  redundancy  for  reliability. 

6.  Using  the  simulation  "fly-off"  results,  attach  rough  cost 
estimates  to  Shuttle  launches,  Upper  Stage  &  EQTV  vehicle  purchases, 
and  payload  delivery  operations.  Compare  coarse  life-cycle  costs 
between  non-reuseabl e  and  reuseable  systems. 

7.  Use  the  overall  optimization  /  methodology  to  investigate 
other  missions  such  as  on-orbit  share  placement  and  recall,  spent 
satellite  retrieval  and  refurbishment,  and  satellite  visitation  for 
refueling  RCS,  battery  /  sensor  replacement,  or  intelligence  gathering 
on  unfriendly  systems. 

These  represent  the  specific  objectives  of  the  research  effort. 
In  the  next  section,  the  approach  chosen  to  meet  these  objectives  and 
justification  for  the  approach  are  presented. 


Approach  and  Justification 

The  approach  and  rationale  for  the  optimization  /  trade-off 
studies  for  different  electric  thruster  technologies  will  be  discussed 
first.  Then  the  approach  and  rationale  for  the  simulation  will  be 
discussed.  The  reason  for  choosing  to  do  both  an 
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is  that  none 


optimization/trade-study  and  a  “fly-off" 
references  which  this  author  has  studied  have  done 
decision-makers  would  like  to  compare  not  only  thruster 


of  the 
both.  Yet, 
technologies , 


but  also  examine  operational  performance  and  costs  over  the  long  term. 


say,  20  years.  That  is,  when  considering  how  to  make  better  use  of 


the  Shuttle  for  delivering  payloads,  and  when  considering  the  benefits 
of  reuseabi 1 i ty ,  one  would  like  to  first  compare  alternatives  and  then 
take  an  optimized  system  and  “fly"  it,  even  if  only  in  a  simulation, 
against  the  baseline  upper  stages. 

Electric  propulsion,  as  mentioned  before,  offers  much  more 
efficient  use  of  propellant  than  chemical  combustion  or  nuclear 
heating/i sen  tropic  expansion.  It  imparts  significantly  more  energy  to 
each  particle  of  propellant  and  greatly  reduces  the  propellant  mass 
needed  to  accomplish  a  given  mission.  This  higher  exhaust  velocity 
provides  mission  capabilities  not  possible  before,  such  as  retrieval 
of  satellites  from  geosynchronous  orbit,  or  delivering  several 
satellites  to  destination  orbits  without  refueling.  iiie  ^ade-off  is 
in  power  supply  size/mass  and  in  very  low  thrust  (hence  long  transfer 
times  between  orbits).  The  long  transfer  time  could  pernaps  be 
tolerated  by  planning  the  launch  date  earlier  and/or  by  hardening  the 
satellite  and  vehicle  against  Van  Allen  radiation.  But  in  all  cases 
for  EP,  it  is  possible  and  desirable  to  shorten  the  transfer  time 
through  choice  of  thruster  and  through  vehicle  optimization. 

The  approach  chosen  for  the  optimization  can  be  summed  up  very 
broadly  in  two  words,  "minimize  mass."  In  doing  this,  cost  generally 
decreases  both  for  earth  launch  (f/kg  to  LEO),  and  also  in  materials 


and  manufacturing  (many  references  express  costs  in  terms  of  S/'kg  of 
the  finished  vehicle).  If  the  total  vehicle  mass  is  fixed,  then 
finding  the  minimum  for  propellant  and  power  supply  mass  means  more  of 
the  total  can  be  payload.  If  all  other  things  were  constant,  transfer 
time  would  also  be  least  for  minimum  mass.  But  the  salient  mass 
trade-off  (between  propellant  consumed  and  power  supply  mass)  is  a 
function  of  the  exhaust  velocity  or  specific  impulse  of  the  thruster 
which  is  governed  by  design  and  by  input  electrical  power.  As  ISP 
increases,  total  propellant  used  for  a  given  mission  decreases,  power 
supply  mass  increases,  and  thrust  increases,  all  of  which  means  that 
transfer  time  continues  to  decrease  with  increasing  input  power  (and 
resulting  thrust  increase) (39,48,41)  .  This  relationship  of  input 
power  to  thrust  and  exhaust  velocity  is  true  for  electron  bombardment 
ion  thrusters  but  not  necessarily  for  other  types  of  EP.  For 
instance,  a  Martin  Marietta  study  (44)  indicates  that  increased  input 
power  is  accompanied  by  decreased  ISP  and  by  increased  thrust  for  the 
self-field  magnetopl asmadynamic  (MPD)  thruster.  As  long  as  a 
mathematical  relationship  can  be  drawn,  though,  between  relevant 
parameters,  the  potential  exists  for  optimization. 

Because  many  of  the  mathematical  relationships  are  non-linear,  the 
optimization  technique  chosen  is  the  Sequential  Unconstrained 
Minimization  Technique  (SUMT)  for  nonlinear  programming.  It  was 
recommended  by  faculty  having  extensive  experience  with  its 
flexibility  in  handling  a  wide  range  of  problems.  At  first  the 
approach  was  going  to  follow  Stuhlinger's  optimization  equations  (52) 
and  maximize  payload  ratio,  maximize  terminal  velocity,  and  minimize 
transfer  time  simultaneously.  But  after  initial  formulation. 


de-bugging  and  examination  of  the  results,  the  output  parameters  were 


rV 


v"' 

'%<• 


m 

m 


fc* 


r* 

N 

tv 

[v  V 


I 


cS 

K: 


fat 


Q 


v:- 


for  a  purely  mathematical  model  and  did  not  represent  any  existing 
thruster.  To  rectify  this  would  require  extensive  modifications  to 
the  equations  by  adding  thruster  efficiencies  and  actual  thruster 
relationships  of  input  power,  exhaust  velocity  and  thrust. 

The  approach  finally  decided  upon  was  inspired  by  Jahn's  text 
(24).  He  suggested  that  as  ISP  increases,  power  supply  mass  increases 
and  propellant  mass  decreases  for  a  given  specific  power  of  the  power 
supply.  Rather  than  aiming  for  the  highest  ISP  in  a  design, 
theoret ical 1 y  an  optimum  ISP  should  exist  for  which  the  major 
components  of  mass  which  vary  (propellant  and  power  supply)  would  sum 
to  a  minimum.  However,  no  method  of  finding  this  optimum  was  given  in 
the  text.  This  relationship  is  also  mentioned  in  the  NASA -us.  is 
1 i terature  (42) . 

With  the  idea  of  minimizing  mass,  it  was  then  necessary  to  find 
the  relationship  of  input  power  to  ISP  and  thrust.  Experiments 
performed  at  NASA-Lewis  Research  Center  (LeRC)  provided  data  points 
from  actual  operating  electron  bombardment  ion  thrusters  of  different 
configurations  (39,40,41,45).  These  thrusters  are  derivatives  of  the 
highly  developed  baseline  30-cm  Kaufman  ion  thruster,  J-series.  This 
thruster  represents  the  closest  to  flight-ready  of  any  primary 
propulsion  EP  design. 

The  data  points  were  first  plotted  and  it  was  noted  that  the 
curves  were  somewhat  linear.  A  linear  curve-fit  was  applied  and 
relationships  for  input  power  vs  ISP  (k  exhaust  velocity)  and  input 
power  vs  thrust  were  developed.  The  first  became  an  input  equation 
for  SUMT .  The  important  added  benefit  of  these  relationships  is  that 
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ail  the  eff iciencies  (mass  utilization,  electric  to  beam  power,  etc) 
are  included  /  accounted  -tor,  since  the  data  is  measured  or  derived 
from  actual  hardware  performance.  With  other  standard  relationships 
from  the  rocket  equation  and  a  relationship  for  specific  power 
(characterizing  the  solar  arrays)  included,  the  SUMT  program  was 
formulated  and  de-bugged.  Further  details  on  equations,  outputs,  and 
the  program  itself  are  found  in  chapter  4. 

The  approach  to  the  simulation  stemmed  from  an  earlier  study  by 
this  author  demonstrating  the  benefits  of  reuseability  of  upper  stages 
/  tugs  in  a  QGERT  simulation  model.  In  this  thesis,  similar  models 
are  developed  for  each  of  the  candidate  vehicles:  Electric  Orbit 
Transfer  Vehicle  (EOTV) ,  Inertial  Upper  Stage  (IUS),  Centaur-G 
(CNTAR),  and  a  representative  Reuseable  Bi-Propellant  Vehicle  (RBPV) . 
The  EOTV  is  optimized  in  StfTT  for  the  payload  or  mission  model  and  is 
reuseable.  The  IUS  and  CNTAR  represent  our  current  upper  stage 
capability  for  heavier  satellites.  For  lighter  satellites,  Payload 
Assist  Modules  (PAMs)  and  kick  motors  are  used.  The  RBPV  is  a  Centaur 
derivative  using  LH,  LOX  in  the  RL-19  engine.  Its  general 
characteristics  are  combined  from  a  Boeing  study  (15)  and  a  Systems 
Engineering  study  performed  at  AFIT  (4?) .  The  proposed  RBPV  does  not, 
however,  include  a  ballute  or  use  aero-braking  as  in  the  case  of  the 
the  Boeing  OTV.  It  is  assumed  to  be  a  space-based  vehicle  requiring  a 
special  refueling  pallet  and  astronaut  /  specialist  team  for  the 
refueling  operation. 

Although  the  SLAM  simulation  program  was  first  considered  due  to 
its  flexibility,  it  was  determined  that  QGERT  (37)  would  be  sufficient 


for  this  study.  Also,  this  author  was  more  familiar  with  Q6ERT  and 


had  used  it  for  a  similar  modeling  problem.  QGERT  allows  modeling  of 
performance  over  selected  periods  (eg.  18,  28,  38  years)  and  allows 
the  modeling  of  uncertainty  to  reflect  real  world  operations. 
Uncertainty  can  be  included  in  reliability  factors,  in  Shuttle  launch 
schedule,  and  any  other  contingencies  desired  by  the  decision-maker. 
Modeling  these  uncertainties  as  well  as  variations  in  schedules, 
changes  to  payload  manifest,  GTV  refurbishment,  etc.,  would  normally 
be  difficult  to  accomplish  analytically. 

The  simulation  “fly-offs"  between  each  vehicle  are  run  over  a  28 
year  period.  The  number  of  shuttle  launches  required  and  the  mass 
delivered  to  orbit  are  determined.  For  the  EGTV,  the  average  transfer 
time  is  obtained.  Inputs  to  the  program  are  the  mission  payload  mass, 
velocity  increment  representing  the  final  orbit,  and  for  the  EGTV ,  the 
equations  characterizing  the  thruster  and  vehicle.  The  outputs  of 
Shuttle  launches,  numbers  of  vehicles  (each  type)  needed,  and  other 
considerations  are  then  used  as  a  basis  for  assigning  coarse  costs  and 
comparing  alternative  operations  during  the  "fly-off." 


Methodol oov 

This  section  will  deal  first  with  a  broad  picture  of  the 
methodology  and  then  show  a  more  detailed  view  of  the  tasks, 
assumptions,  and  information,  fl ow  that  are  involved  in  doing  analysis 
with  this  methodol ogy. 


It  can  be  seen  in  Figure  1  that  the  mission  determines  the  needs 


and  provides  inputs  tor  the  3UMT  and  QGERT  programs.  Each  mission  is 
characterized  by  three  major  parameters:  mass,  orbit  altitude,  and 
orbit  inclination.  The  desired  orbit  altitude  and  inclination  are 
used  to  determine  the  required  velocity  increment,  Av,  through  use  of 
the  A1 f ano-Wiesel  curves  (1).  These  curves  give  the  minimum  energy 
transfer  in  the  case  of  low  thrust  (ie.  -  the  EQTV)  .  Use  of  these 
curves  is  also  covered  in  Appendix  II.  The  Av  -for  the  other 
vehicles,  IUS,  CNTAR,  RBPV,  is  a  straight-forward  Hohmann  transfer  and 
because  these  transfers  take  place  in  less  than  a  day,  no  detailed 
calculations  are  required  by  this  methodology.  Therefore,  for  QGERT 
"fly-off"  purposes,  chemical  systems  are  assumed  to  make  the  transfer 
in  0.5  day. 

The  methodology  is  set  up  so  that  the  user  may  specify  a  mission 
or  a  mission  set  and  run  the  algorithms^  with  this  input.  The 
algorithm  is  designed  primarily  for  a  deployment  scenario  in  which 
payloads  are  transported  from  LEO  to  some  higher  orbit  and 
inclination.  Other  scenarios  can  be  envisioned  and  these  are  treated 
in  Appendix  V.  For  discussion  purposes,  the  NavStar  GPS  has  been 
selected  as  a  representative  mission  payload  for  use  through  the 
entire  methodology.  Chapters  4  and  5  will  discuss  the  specific  inputs 
and  outputs  in  some  detail,  but  the  emphasis  is  on  the  overall  oicture 
in  this  chapter. 

Referring  to  Figure  1  again,  the  lower  left  block  depicts  the 
optimization  of  an  electric  propulsion  Orbit  Transfer  Vehicle  (EOTV) . 
Components  chosen  for  the  propulsion  and  Dower  systems  come  mostly 
from  the  30-cm  Ion  Thruster  Design  Handbook  (12)  plus  a  compilation  of 
data  from  numerous  other  sources  in  the  bibliography.  The  user  mav 


select  candidate  EP  propulsion  systems  as  desired,  as  long  as  the 
relationship  between  input  power  and  ISP/thrust  can  be  ex Dressed 
mathematically.  It  not,  the  user  must  modity  the  SUMT  program  incut 
more  extensively.  For  the  GPS  example,  five  thrusters  compatible  with 
the  Design  Manual  were  chosen  and  optimized  tor  1 ,2, 3, 4, 5, and  6,  GPS 
satellites  aboard.  Based  on  the  minimum  mass  ot  the  system  Determined 
by  the  SUMT  optimization  and  on  secondary  calculations  ot  the  transter 
time,  an  acceptable  combination  ot  minimum  mass  and  minimum  transter 
time  are  chosen  as  the  desired  optimum  EQTV.  Details  will  be  covered 
in  Chapter  4. 

Having  chosen  an  ootimal  EOTV,  this  vehicle  can  then  be  examined 
tor  long  term  pertormance  in  a  "tly-ott3  against  our  present  upper 
stages  and  a  representative  Reuseable  Bi-Propellant  Vehicle  (RBPV)  in 
the  QGERT  simulation.  This  phase  is  shown  in  the  lower  right-hand 
block  in  Figure  1.  Mission  characteristics  ot  Av  required,  Shuttle 
integration  tactors,  etc.,  and  vehicle  characteristics  tor  EOTV  are 
inputs  tor  the  QGERT  models.  One  QGERT  model  exists  tor  each  ot  the 
upper  stages  and  OTVs.  The  user  may  decide  the  appropriate  time 
period  to  run  the  simulation,  but  20  years  is  used  in  this  thesis. 
While  the  simulation  may  be  run  using  a  single  type  ot  satellite 
payload,  a  mission  set  can  be  used.  Either  may  be  specitied  by  the 
user.  More  details  will  be  covered  in  chapter  5. 

With  the  results  ot  the  simulation,  the  pertormance  ot  each 
vehicle  may  be  examined  and  costs  assigned  to  the  -operations. 
Determining  costs  has  proven  to  be  ditticult  and  elusive  tor  systems 
not  yet  in  existence,  and  even  tor  the  IU3  since  it  has  not  had  as  big 


a  block  buy  as  originally  intended.  However,  coarse  cost  estimates 


can  be  used  -for  comparison. 


Other  results  o-f  interest  to  the 


decision-maker  are  the  nunoer  of  Shuttle  launches  required  in  using 
each  upper  stage  /  G7V  and  the  total  number  o-f  satellites  delivered, 
y-lith  these  results  and  a  knowledge  o-f  the  assumptions  and  inputs 
comprising  the  results,  the  decision-maker  should  be  in  a  better 
position  to  evaluate  the  merit  o-f  operating  with  one  or  a  mixture  o-f 
tnese  vehicl es. 

To  grasp  more  o-f  the  details  in  working  with  this  methodology, 
refer  to  Figure  2  for  the  following  discussion.  The  operational  need 
drives  the  requirement  for  a  satellite  system.  Those  satellites  not 
deployed  in  LEO  must  be  delivered  to  higher  orbit  by  kick  stages  or 
Payload  Assist  Modules  (PAMs)  for  lighter  satellites,  or  IU3, 
CENTAUR-G  for  heavier  ones.  The  decision-maker  then  has  a  need  for  a 
tool  or  method  for  evaluating  these  expendable  stages  against  more 
fuel -efficient ,  reuseable  electric  OTVs  and  against  reuseable  chemical 
OTVs. 

For  purposes  of  the  optimization,  the  user  may  choose  a  specific 
payload  or  a  representative  mission  set.  If  a  specific  payloao  is 
chosen,  the  SUMT  program  will  require  the  mass  of  the  payload  to  be 
inserted  in  the  equation  tableau.  A  comment  card  is  included  by  this 
author  at  the  appropriate  lines  in  the  formulation.  If  a  mission  set 
is  chosen,  three  approaches  exist.  The  first  is  to  run  SUMT  ano 
optimize  an  EQTV  for  every  different  payload  in  the  set.  A  single 
optimum  EQTV  is  then  chosen  to  favor  the  payload  with  the  highest 
frequency  of  launch,  or  a  weighting  technique  can  be  devised  with 


weighting  factors  determined  by  the  decision-maker. 


The  second 
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FIGURE  2.  DETAILED  METHODOLOGY 


approach  is  to  categorize  the  whole  set  by  a  single  reoresen t i t ive 
payload  and  optimize  the  EOTU  for  that  nominal  payload.  The  third 
approach  is  to  group  the  set  into  similar  mass  and  orbit  requirements 
and  optimize  a  fleet  of  EOTVs,  one  for  eacn  group. 

The  orbit  altitude  and  inclination  for  the  payload  must  be  used  in 
the  A) f ano-Wiesei  transfer  to  obtain  the  Av  required.  This  Av,  in 
KM/sec,  is  an  input  for  SUMT  and  for  subsequent  transfer  time 
calculations.  It  appears  in  two  equations  in  the  SUMT  formulation  and 
is  readily  edited.  With  Av,  payload  mass,  and  general  EOTV 
characteristics  as  inputs,  the  SUMT  formulation  allows  the  user  to 
choose  an  appropriate  thruster  technology  and  optimize  the  EOTV.  For 
^nitial  runs,  the  general  EOTV  parameters  were  masses  from  the  Design 
Manual  (12).  The  selected  thrusters  were  all  derivatives  of  the  30-cm 
electron-bombardment  Kaufman  thruster.  The  user  may  optimize  given 
payloads  with  these  parameters  or  select  other  EOTV  configurations  and 
thrusters.  SUMT  is  primarily  minimizing  the  combined  masses  of  the 
required  propellant  and  required  power  supply.  In  the  ootimization , 
the  parameters  that  are  being  determined  by  SUMT  are  the  throttling 
level  for  the  set  of  thrusters  (input  power  for  given  exnaust 
velocity,  thrust,  and  mass  flow)  and  the  resulting  trade-off  between 
fuel  mass  consumed  and  power  supply  mass  required  to  do  the  mission. 
The  mission  is  two-way,  deploy  and  return  with  most  of  the  fuel 
consumed.  Note  that  electric  propulsion  systems  have  the  potential 
for  accomplishing  several  "ou t-and-bacKs“  without  refueling,  whereas 
chemical  propulsion  systems  are  hard-pressed  to  make  one 
"out-and-back . 3  The  optimization  results  provide  the  required  input 
power  to  the  thruster  in  order  to  simultaneously  minimize  the  fuel 
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mass  and  power  supply  mass. 


Since  the  structure. 


housekeeping,  and 


payload  masses  have  been  specified,  the  major  drivers  are  prooellant 
and  power  supply.  The  optimum  EQTV  has  been  -found  by  3UMT  as  iong  as 
the  associated  transfer  time  is  acceptable. 

The  equation  to  determine  round-trip  transfer  time  for  the  optimum 
EOTV  deploying  its  associated  payload  is  given  by: 


rp  _  (Mpi  A1  pz)  3c  (isp) 

H  ("Hi) 


i a  round-trip  transfer  in  seconds,  Mpj  is  mass  of  propellant  out 
to  orbit,  M  p2  is  mass  of  propellant  return,  N  is  number  of  thrusters, 
and  Th  is  thrust  per  thruster.  All  units  are  MKS.  This  is  a  standard 
relationship  based  on  a  rearrangement  of:  thrust  =  mass  flow  x 
exhaust  velocity.  Masses  for  the  equation  are  obtained  from 
derivatives  of  the  rocket  equation,  eq.  1.  If  the  transfer  time  is 
not  acceptable,  more  thrusters  can  be  added  and  the  EOTV  optimized 
again,  or  more  input  power  can  be  supplied  (away  from  the  minimum 
mass)  until  the  transfer  time  is  within  an  acceptable  range.  CIt 
should  be  mentioned  here  that  a  program  was  formulated  which  minimized 
both  mass  and  weighted  transfer  time.  However,  the  solution  from  SUMT 
would  not  converge  within  central  processor  time  limits.  That  effort 
is  recommended  for  further  study.] 

Once  the  transfer  time  is  acceptable,  the  methodology  flows  to  the 
"fly-off“  simulation.  The  characteristics  of  the  optimal  EOTO,  the 
payload  mass(es)  and  Av(s)  ,  and  any  special  mission  constraints,  such 
as  Shuttle  compatibility  /  integration,  transfer  time  restrictions, 
etc.  are  input  to  the  QGERT  models.  One  model  for  each  upper  stage  / 
transfer  vehicle  was  developed.  This  allows  a  clearer  comparison 
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between  them  as  each  delivers  the  same  representative  payload  or 
mission  model .  The  characteristics  tor  IU3,  CNTAR,  and  RBPV  whicn 
were  previously  determined  are  incorporated  in  each  model . 

Each  of  the  tour  vehicle  QGERT  models  is  directly  intluenced  bv 
the  Shuttle  launch  rate.  Because  the  RBPU  is  requiring  extra  Shuttle 
tlights  to  bring  up  tuel  and  because  valuable  payload  bay  space  is 
taken  by  the  upper  stage  in  the  case  ot  expendable  systems,  there  is 
potential  tor  reducing  the  Shuttle  launch  rate  with  an  tOTV. 
Requisite  tor  each  simulation  model  is  the  turn-around  tor  the  Shuttle 
—  how  many  trips  to  LEO  are  possible  in  a  year?  Given  the  ditticult 
tasks  ot  returbishment ,  ET  and  SRB  mating,  payload  integration, 
retueling,  and  system  checkouts  tor  Shuttle,  the  turn-around  with  3  or 
4  Orbiters  will  probably  not  be  under  20  days  very  otten.  So,  a 
realistic  figure  of  13  launches  per  year  is  chosen  for  late  1939s, 
early  1999s,  even  if  it  is  a  bit  conservative  compared  to  many 
literature  sources.  Several  studies  assume  very  high  launch  rates 
commensurate  with  greatly  increased  space  activity  in  the  1999's  and 
beyond.  But  Federal  and  industry  budget  constraints  are  likely  to 
persist  in  the  next  29  years,  and  that  will  undoubtedly  curtail  both 
customers  and  operators  of  the  Shuttle.  The  lower  assumed  launch  rate 
may  tend  to  favor  the  expendable  systems,  and  therefore  the  study 
results  should  be  conservative  in  that  regard,  too.  Examining  total 
number  of  satellites  delivered  in  the  20  year  (selectable  by  user) 
period  and  comparing  to  the  number  required  by  the  mission  model  can 
suggest  what  reduction  in  Shuttle  launch  rate,  if  any,  is  possible. 

QGERT  readily  allows  for  introduction  of  uncertainty.  While  a 
decision-maker  may  want  to  mode!  uncertainty  throughout  the  "fly-off", 
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this  author  chose  3  areas  that  in  reality  should  nave  uncertainty 
associated  with  them.  The  -first  is  the  input  Shuttle  launch  rate.  A 
no.-~al  distribution  is  chosen  with  standard  deviation  or  +/-  davs. 
An  earlier  modeling  nad  incorporated  an  uneven  interval  between 
launches  with  an  exponential  interarrival  distribution.  But  the 
inherently  large  variance  was  not  needed  in  this  “tly-off." 

The  second  area  o-f  uncertainty  is  -failure  o-f  the  upper  stage  or 
trans-fer  vehicle.  For  the  EGTJ,  a  reliability  o-f  .999  is  assumed 
because  i-f  one  or  even  several  thrusters  -fail,  the  mission  is  not 
aborted.  The  trans-fer  will  simply  take  longer  as  the  remaining 
thrusters  are  gimballed  to  compensate  -for  the  lost  thrusters  and  the 
mission  continued.  At  low  thrust  levels  this  is  possible,  whereas  ail 


chemical 

vehicles  being  compared 

here  would 

be 

total  1 y 

1  ost 

upon 

t'ai  1  ure 

of  one  engine  (the  on  1 v 

engine 

in 

the 

case  of 

IUS 

kick 

stages)  . 

Also,  the  EOTV  has  at 

least  two 

PPUs 

per 

BIMOB 

unit, 

two 

inter-face  modules  per  vehicle  and  at  least  two  solar  panels.  It  is 
assumed  that  the  malfunctioning  unit  could  then  be  replaced  upon 
return  to  Shuttle  orbit  and  arrival  of  the  next  Orbiter  with  parts. 
Assumptions  for  other  model  reliabilities  are:  .965  for  IUS,  .935  tor 
Centaur  (weighted  by  previous  performance)  and  .975  for  the  RBPU.  The 
decision-maker  is  free  to  make  alternate  assumptions  for 
rel iabi I i ties. 

The  third  area  of  uncertainty  is  in  Shuttle  payload  manifest, 
since  changes  will  undoubtedly  occur  to  the  schedule.  Also,  some 
payloads  will  not  require  transfer  to  higher  orbits.  To  account  for 
this,  a  node  is  included  where  65/i  of  the  time  the  arriving  pavload  is 
not  transferred  by  any  of  the  four  vehicles.  Reasons  for  this  might 


include,  Spacelab  mission,  PAM  deployed  payload,  LEO  experiment 
package,  and  LEO  satellites.  More  will  be  discussed  in  Chapter  5.  It 
is  important  to  note  that  although  uncertainty  is  included,  each  o f 
the  -four  QGERT  models  will  be  subjected  to  the  same  random  number- 
stream  so  that  variance  between  compared  runs  is  controlled.  That  is, 
uncertainty  affects  outcomes  0+  the  -four  models,  but  is  applied 
consistently  to  each  model. 

With  all  of  these  results  of  reduced  Shuttle  launch  rates,  OTU 
fleet  sizes,  number  of  satellites  delivered  to  final  orbits,  and 
effects  of  uncertainties,  the  methodology  proceeds  to  the 
determination  of  coarse  Life  Cycle  Costs  (LCC)  for  acquisition  and 
operation  of  each  alternative  vehicle.  More  will  be  discussed  as  to 
assignment  of  costs  in  Chapter  4.  With  the  performance  of  each  OTU  / 
Stage  and  a  rough  comparison  of  LCCs,  the  decision -maker  now  has  a 
tool  with  certain  flexibilities  that  allow  comparison  and  evaluation. 
With  more  information,  the  decision-maker  can  assess  a  decision  to 
continue  present  operations  with  IUS  and  later  the  Centaur-G  or  to 


acquire  either  electric  or  chemical  reuseable  transfer  vehicles. 


CHAPTER  IU.  Optimization;  Using  SUMT 


SUMT  Non-Linear  Program 

The  Sequential  Unconstrained  Minimization  Technique  (SUMT)  program, 
has  been  chosen  to  optimize  EGTVs  -for  given  payloads.  SUMT  has  tne 
flexibility  to  use  more  than  one  available  techmaue  to  find  tne 
minimum  of  a  multivariable,  nonlinear  function  subject  to  nonlinear 
inequality  and  equality  constraints. 

The  genera]  mathematical  programming  formulation  of  the  oroblem  is: 

Minimize:  F<Xi>  Xi,...,  Xw) 

Subject  to:  6k(Xlf  X2,...,  Xri)  >8,  k=  1,  2,  M 

Hk(X1,  X2 ,  .  . .  ,  XH>  *0,  k=  M+l,  M  +  2,  ....  M+M2 

The  procedure  for  solving  this  non-linear  programming  (NLP)  problem  was 
developed  by  Fiacco  and  McCormick  as  detailed  in  references  20.  33  and 
57.  The  inequality  and  equality  constraints  are  attached  to  the 
original  objective  function  to  make  use  of  penalty  function  techniques. 
The  resulting  unconstrained  function  is  minimized  using  one  of  several 
appropriate  mu  1 tivariabl e,  unconstrained  techniques. 

The  following  quote  and  diagram  (Figure  3)  from  reference  29 

describe  the  procedure  followed  by  the  algorithm  in  finding  the 

minimum. 

"The  algorithm  proceeds  as  follows: 

1)  A  modified  objective  function  is  formulated  consisting  of  the 
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original  function  and  penalty  functions  with  the  form: 


M  M-MZ  _ 

P.P-r  S  Aifi  +  E  H C/r 
k-1  kxM+i 

where  r  is  a  positive  constant.  As  tne  algorithm  progresses,  r  is 
reevaluated  to  form  a  monotonical 1 y  decreasing  sequence  r1  >r2  >...>0. 
As  r  becomes  small,  under  suitable  conditions  P  approaches  F  and  the 
problem  is  solved. 

2)  Select  a  starting  point  (feasible  or  nonfeasible)  and  an  initial 
value  for  r. 

3)  Determine  the  minimum  of  the  modified  objective  function  for  the 
current  value  of  r  using  an  appropriate  technique  (several  options 
ava i 1 ab 1 e) . 

4)  Estimate  the  optimal  solution  using  extrapolation  formulas  (33) . 

5)  Select  a  new  value  for  r  and  repeat  the  procedure  until  the 
convergence  criterion  is  satisfied."  (29) 

The  developers  of  the  code  consider  it  to  be  a  research  tool  ratner 
than  a  production  code  for  NLP  (33).  This  is  because  it  allows 
experimentation  with  various  techniques  for  solving  NLPs,  since  no 
existing  technique  handles  all  types  of  problem  formulations.  It  is 
this  flexibility  that  has  made  it  a  choice  for  optimizing  EGTVs. 


Preliminary  Mission  Model 

It  should  be  made  clear  that  the  user  of  this  methodology  is  able 
to  choose  a  single  payload  or  several  payloads  comprising  a  mission 


model  and  then  have  the  EOTV  optimized  for  that  payload  or  model. 


For  purposes  of  demonstrating  and  verifying  this  portion  of  tne 
methodology,  the  NavStar  GFS  satellite  constellation  has  been  selected 
as  a  mission  model.  This  allows  preliminary  results  to  be  obtained  tor 
thesis  sponsors  at  USAF  Space  Division.  Also,  other  payloads  may  be 
represented  by  the  mass  categories  arising  from  deploying  1  -  6  GPS 
satellites  at  the  same  time  with  one  OTU.  For  instance,  one  FLTSATCOM 
satellite  would  be  similar  in  weight  to  2  next-generation  GPS 
satel 1 i tes. 

The  -following  table  shows  the  mission  model  in  terms  o-f  number  o-f 
GPS  satellites  (3AT3) ,  mass  in  KG,  weight  in  LB,  and  Av  -for  LEG  to  GEO 
orbit  change.  Current  GPS  satellites  are  in  sub-synchronous  orbits, 
but  GEO  is  a  possibility  in  the  -future.  So,  this  "worst  case"  scenario 
of  deploying  from  a  209  km  orbit  to  GEO  has  been  used  during  initial 
optimization  runs. 


xmmmmmxmxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxmmxxxxx*s*xs**s* 

TABLE  1 


#  of  SATS 

MASS  <KG) 

WT  (LB) 

Ay  to  GEO  (KM/SEC) 

1 

993 

2092 

5.3332 

2 

1816 

4904 

5.3332 

3 

2724 

6905 

5.3332 

4 

3632 

3907 

5.3332 

5 

4549 

10009 

5.3332 

6 

5443 

12911 

5.3332 

The  Av  in 

Table  1  was 

obtained  from 

the  Al f ano-Wi esel  Optima 

Many-Orbit  Transfer  Curves  for  I  ow  thrust.  The  weight  figures  resulted 
from  rounding  the  mass  of  the  2900  lb  ol anned-growth  GPS  to  993  Kg. 


Multiples  of  ?9S  were  used  for  the  runs  and  comprise  the  KG 


co i umn  . 


t.-'' 


When  converted  back  to  lbs.  -for  comparison,  the  weights  are  not  'nice 
and  round."  MKS  units  are  used  throughout  this  methodology  tor  all 
calculations.  "PI anned-growth “  refers  to  the  tact  that  other  sensors 
and  hardening  are  planned  tor  tuture  GPS  satellites.  These  satellites 
should  be  ready  tor  deployment  in  the  time  frame  that  EOTVs  coulo  reach 
Initial  Operational  Capability  (IOC). 


Characterizing  the  EQTU 

In  order  to  make  initial  runs,  assumptions  have  been  made  regarding 
Electric  Propulsion  Orbit  Transfer  Uehicle  masses  and  these  assumptions 
will  be  discussed  first.  Next,  the  method  for  characterizing  the 
thrusters  for  the  SUMT  program  will  be  discussed.  For  more  information 
on  the  proposed  vehicle  configuration,  the  reader  is  referred  to 
Appendix  I. 

The  30-cm  Design  Manual  (12)  contains  tables  of  masses  for  the 
components  of  the  thruster  subsystem  called  B I MOD .  As  the  name  lmolies 
this  BIMOD  unit  is  modular  and  includes:  two  thrusters  and  their 
associated  Power  Processing  Units  (PPUs) ,  thermal  control  for  PPUs,  and 
thruster  gimbals.  When  fastened  together,  4  BIMOD  units  containing  3 
thrusters  form  the  basic  propulsion  configuration  assumed  for  the 
preliminary  optimization  runs.  On  top  of  the  4  BIMOD  units  is  an 
interface  modul e  which  contains  propellant  tanks,  gimbal  electronics. 


thruster  controller,  power  distribution,  truss  structure,  harness,  and 
more  thermal  control.  The  given  mass  of  each  BIMOD  unit  is  137.1  Kg 
and  thus  4  units  total  548.4  Kg.  The  single  interface  module  for  these 
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4  BIMGDs  is  153.7  Kg,  less  propellant.  The  total  propulsion  suosvstem 
is  797.1  Kg.  Telemetry,  guidance  and  control,  and  other  avionics  pius 
Control  Moment  Gyros  (CMGs)  in  lieu  of  a  Reaction  Control  System  (RCS) , 
solar  panel  array  steering,  autonomous  computing,  and  pavload 
inter-faces  are  assumed  to  round  the  vehicle  mass  to  1000  Kg.,  total. 
Foliow-on  runs  will  include  -further  assumptions  as  to  vehicle 
characteristics.  Users  are,  of  course,  free  to  make  their  own 
assumtions  /  calculations  of  the  EGTV  masses  as  new  data  is  avail aole. 

With  an  assumed  vehicle  mass  of  1000  Kg.,  the  next  task  is  to 
determine  a  way  to  characterise  thruster  performance  in  the  SUM! 
program.  This  has  been  accomplished  by  noting  that  experimental  data 
on  ion  thrusters  (39,40,41,45)  shows  an  aoprox imatel y  linear 
relationship  between  thruster  input  power  ana  performance 
characteristics  of  ISP  (  &  exhaust  velocity  )  and  thrust.  Figure  4  is 
a  graph  showing  ISP  and  thrust  plotted  against  input  power  for  a 
Ring-Cusp  3-grid  ion  thruster.  As  input  power  is  increased,  both 
ordinate  values  increase.  A  simple,  linear  curve  fit  was  performed  for 
each  of  5  electron-bombardment  ion  thrusters.  Examples  of  the 
experimental  data  table  and  ISP  and  thrust  curve  fits  are  in  Table  2 
and  Figure  5,  respectively.  These  data  points  and  the  curve  fit  are 
for  the  Ring-Cusp  3-grid  ion  thruster.  Data  and  curve  fits  for  the 
other  4  thrusters  are  found  in  Appendix  I.  In  follow-on  studies  it  is 
hoped  to  include  other  electric  propulsion  technologies  such  as 
self-field  magnetopl asmadynamic,  pul sed-inductive,  rail-gun,  and  other 
promising  designs. 


FIGURE  4. 


Tnpiat  Power,  vs  xsp >  thrust 
RING)  -  cusp  3 -Gem  xeajon  thruster. 


PROJECTED  RING-CUSP  THRUSTER  PERFORMANCE  WITH  INERT  GAS  PROPELLANTS  USING  TWO  TYPES  Of  ION  OPTICS 
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TABLE  Z.  Ring  -  cusp  data .  from  sjas a  (is). 


TABLE  -  Tsp 


IN  RUT  P **■>  Wj  M 

PATH  ISP 

Cutvt  Fir 

0/FF£P£A/C£ 

EM 

/.&  KW 

1580 

/ 58 8.  8 S' 

s.fsox 

.00*/ 

2.0C 

20 SO 

1033. 29 

C.7 08 

.00? 

2. 32 

2Z70 

22f0.?57 

10.7S7Z 

.00s 

JOsr  STATS  i  M£aH  Nt  J.W  °Mx~’  *a 2,03 7  5*  =  0,37 

M£AU  XsPi  rtfc'7  <tNy  r  287.  ms  Sy  * 352.  #7S 

EQtO  TLst  =  fs/.  7 SIC  (w)  +  72.  Co  S3 

SUMT  FoZHULATJotJ:  3/  VAL  =  Xfi)  -  957,  7*  ft  (w)  -  7Z .  C0S3 


TABLE  -  THRUST 

INPUT  PwR,  (tw) 

77/pusr(v) 

cutve  Frr 

J>IFP£<£AICE 

% 

AS?  XL/ 

.09*/ 

.09*/ 139 

.00013 

.00/ 

2.0C 

.121 

.IZ0C09 

000311 

,003 

2.32 

.13S 

.13525/8 

.0002S/  8 

.002 

THRUST  stats  :  MEAN  W:  1.99  crNX  =  ,  302/037  Sx  =  0.57 

At  fa*  Th:  .  1167  *] it  0/?0/t3  Sy  - ,  ozo8io7 

£M:  Thrust  -  .0503/8  (W)  -t.00*/S93 

SUMT  Po/tMULATioHi  3*1  VAL  =  X  fa)  -.OSC319  (w) ~ .009S93 


FIGURE  5. 


THRUST J  ISP  I/S  U/  CAUulat/ohs 
RtHQ  -  CUSP  3-<5*fJ>  XE  THRusTfR. 


Several  of  the  available  options  in  the  SUMT  i ibrary  were  tried 
against  four  of  the  types  of  thrusters  for  EQTVs  carrying  payloads  from 
2080  lbs.  to  12000  lbs.  The  results  of  using  different  options  varied 
from  program  dump  for  lack  of  Central  Processor  (CP)  time  to  sucessful 
runs  with  only  slight  variations  (tenths  of  seconds)  of  CP  time. 
Answers  varied  from  exact  (  ><i  decimal  places  )  to  erroneous.  Thus,  a 
significant  amount  of  time  has  been  devoted  to  finding  a  set  of  options 
that  gives  meaningful  answers  within  reasonable  CP  time  limits. 

Option  keys  allow  the  SUMT  user  to  input  different  convergence 
criteria,  printout  options,  and  problem  linearity,  as  well  as  the 
desired  unconstrained  minimization  technique.  The  option  kevs  which 
seemed  best  for  EOTU  formulations  were: 

NT ( 1) =3  r  value  option  set  to  RHOIN 

NT ( 2)  =  1  automatic  inclusion  of  trivial  constraints,  X  2. ^ • 

NT ( 3)  =  1  standard  printout 

NT(4)=1  final  convergence  determined  on  basis  of  current 
subproblem  solution. 

NT (5)  =2  final  convergence  option 

NT(6)=1  no  extrapolation 

NT(7)=1  subproblem  convergence  option 

NT(8)=1  linearity:  at  least  one  nonlinear  constraint 

NEXQP1=1  option  for  checking  derivatives 

NEX0P2=1  unconstrained  minimization  technique  —  in  this  case,  the 
method  chosen  is  the  generalized  Newton-Raphson  method  as 
handle  indefinite  Hessian  matrices. 


modified  to 


In  some  cases  where  convergence  was  requiring  too  much  CP  time  — 
such  as  -for  the  inert  gas  thruster  —  option  NT<4)  was  set  to  2  so  that 
final  convergence  was  determined  on  the  basis  of  first  order  estimates. 
This  required  typically  1/2  the  CP  time  as  before,  but  at  the  expense 
of  accuracy.  Instead  of  6  decimal  places,  this  option  satisfied 
equality  constraints  to  within  only  3  decimal  places  of  the  optimum 
point.  This  is  not  considered  significant  for  the  EQTV  optimization 
problem. 

Formulation  of  the  set  of  equations  for  the  3UMT  algorithm  was 
straight-forward.  Limits  had  to  be  set  for  the  feasible  region  both 
to  insure  realistic  values  and  to  prevent  spending  CP  time  on  too  wide 
a  search.  These  formed  the  inequality  constraints.  The  equality 
constraints  were  based  on  standard  relationships  derived  from  the 
rocket  equation  for  propellant  mass  used  out  to  orbit  (Mp(  )  and  return 
<M  pa  ) .  Also  included  are  the  relationship  for  input  power  vs  thrust 
previously  developed  and  a  standard  relationship  for  the  specific  power 
of  the  solar  power  supply. 

Figure  6  shows  an  initial  formulation  for  the  Ring-Cuso  thruster 
with  3-Grid  ion  optics.  Limits  on  specific  impulse  (XI)  were  50«  to 
10800  sec.  Limits  on  power  supply  (solar  array)  mass  (X2)  were  10  to 
19000  Kg.  Limits  for  both  propellant  masses  (X3  out  to  GEO  and  X5 
return)  were  10  to  10000  Kg.  Finally,  limits  r  input  power  (X4)  to 
the  thruster  were  0.5  to  20  KW.  The  specific  p  er  equality  (X2  vs  X4) 
assumes  no  cabling  resistance  losses  and  assumes  1.0  for  FPU 
efficiency.  Thruster  efficiencies  are  inherent  in  the  thruster 
relationship  of  input  power  <X4)  vs  ISP  (XI). 

Figure  7  shows  a  orogram  printout  of  the  completed  formulation  for 


a  different  thruster  —  the  Ring-Cusp  3-Grid  Xenon  Ion  thruster.  with 
the  exception  of  input  power,  the  constraints  are  the  same.  Line 
number  31  of  the  orogram  input  is  the  net*)  thruster  relationship. 

Figure  3  shows  the  final  results  printout  for  the  formulation  in 
Fig. 5,  and  is  typical  of  the  output  format  for  the  final  solution  to 
the  minimization  problem.  The  final  value  of  F  represents  the 
minimization  of  the  objective  function.  Recall  that  this  function  is 
the  total  mass  of  the  power  supply  <3  S  X2)  -jjs  the  total  propellant 
mass  out  to  GEO  and  back  (X3  +  X5) .  The  X  values  are  just  the  solution 
values  of  the  X  vector  for  the  minimized  function.  Constraint  values 
are  the  equality  and  inequality  constraints  solved  usino  the  final  X 


vector . 


VAR/ABLE  f>EF/rttTJotfsi 

X(l)  —  JTs f  —  Specific  In  pulse.  (pee) 

X  (l)  -  A/jf'  —  A/d-SS  Z?*'*/*  Supp/y  per  lfiusi+r  (kG) 

X(3)  —  Mrt  -  M*SS  *f  fhpdfemt  use/  ,ut  t4 
/(V)  =  u/  =  Input  Power  per  thruster  (bv) 

X(s)  —  fife  =  A/*#  //  Propel t cut  Usuf  fir  Cvturu  trip  (pi) 

FOKMt4L4T/0R  : 

Fcsjj  M/M/AUZE:  F  -  2  *  X(f)  +  x(i)+  X(s) 

Coa/skA/uts  suss  To:  XO)  -Soo 

XO)  *10,000 
Xfr)  *■  10 
X( z)  *  IO.000 
X(3)  -  10 
X  (3)  —  10,000 
X(1)Z0.1 
X(i)<20 
X(5)  *  10 
x(s)  —  i^OOO 

YO)  -  9SA  7X1  L  XX  Cl)  +  72,  C>OSS 
x(z)  =  X(4)/osz 
X(?)  =  {lOOO-i  s*xc4)* 

(exp  (s.9W/ooWXX6j)-i) 
X(S)  =  (.2000  +  X*X(2))+X$)x 

(exp  (s,tm/ooW*XQ))-i) 

SUMT  XNPUT  EQUAT/OtiS: 

08T.  fed  1  YAL  =  f*X(*)  +X($  tk(s) 

EXAMPLE  CoxaiXAft/r:  YAL  =  X(l)  -500,0 


FIGURE  6. 


Problem  Form  flat /oa/  fog.  sumt~ 
RifJG  -  CUSP  3  -  <SAiL>  XEW»*J  THAustEG 


FIGURE  7 


3UMT  Equation  Input  -  Ring  Cusp,  1 0 O O  Kg 


»r  Prooraro  Main; 

=C  RESTRAINT  PORTION 

*  SUBROUTINE  RESTNT  (IN, UAL) 

*  COMMON/ SHARE/X(  100)  ,OEL<  100)  ,A(  188 , 190)  ,N(5)  , 
+M,MN,NP1,NM1 

IF  (X<  1)  .LE.  1.0)  X(l)=500.0 
IF  (IN')  19,10,20 
10  UAL=  3*X<2)  +  X(3)  +  X(5) 

=  RETURN 

*  29  GO  TO  (21 ,22 ,23,24 ,25,26,27,28,29 ,30 ,31 ,32 ,33,34) , IN 

21  UAL  =  X(  1)  -500.0 
RETURN 

22  UAL  =  10000.0-  X( 1) 

RETURN 

23  UAL  =  X(2)  -  10.0 
RETURN 

24  UAL  =  10900.0  -  X(2) 

RETURN 

25  UAL  =  X(3)  -  10.0 
RETURN 

26  UAL  =  10800.0  -  X(3) 

RETURN 

=  27  UAL  =  X(4)  -  0.5 

RETURN 

23  UAL  =  20.0  -  X( 4) 

RETURN 

29  UAL  =  X(5)  -  10.0 
RETURN 


560=  30  UAL  =  10080.0  -  X<5) 

578=  RETURN 

530=  31  UAL  =  X(l)  -  95 1 . 7S96£X( 4)  -  72.6053 

590=  RETURN 

609=  32  UAL  =  X( 2)  -  (X(4)/(0 .052) ) 

610=  RETURN 

620=  33  UAL  =  X(5)  -  (1900  +  3*X<2) ) * ( EXP ( 5 . 3332/ (8 .0093 1*X(  1) > ) - 1) 

630=  RETURN 

640=C  PAYLOAD  MASS  IS:  1000KG 

650=  34  UAL=X(3)-(2000.0+8*X(2)+X(5))£(EXP(5. 3382/(0. 00981*X< 1)))-1 

660=  RETURN 

670=  OvID 

630=*EOR 

690=  $DATA  N=5 ,M= 10 ,MZ=4 , 

700=  X= 1976. 2, 33. 462, 972, 437, 2. 0,459. 533, 

710=  NT( 2)  =  1  ,NT(5) =2,NEX0P2=  1  *END 


V*  - 


a 


IsSfl 


FIGURE  3 

3UMT  Final  Output  Values  -  Ring  Cusp,  1000  Kg 


After  the  Last  Iteration: 

3950= 

xx*xx*xxxx*x**xxxx 

3930=  FINAL  VALUE  OF  F  =  1 .331 1700  1E+03 

3970= 

3930= 

3990=  FINAL  X  VALUES 
9000= 


9010=  X<  1)  = 
4.780491 19E+02 

3.37917993E+03 

X( 

2) 

s 

7.391 12323E+01 

X  (  3)  = 

9020=  X<  4)  =  3.99938SS9E+00 

9030= 

9040=  FINAL  CONSTRAINT  VALUES 
9050= 

X< 

5) 

2 . 3733033 1E+02 

X( 

9030=  G(  1)  = 
3.391 12823E+01 

3.379 17993E+03 

G( 

2) 

= 

3.  120S2004E+03 

G(  3)  = 

9070=  G(  4)  = 
9.52195038E+03 

9.92308372E+03 

G< 

5) 

= 

4.33849 1 19E+02 

Gt  3)  = 

9030=  G<  7)  = 
2.5783033 1E+02 

3.49933339E+00 

G( 

8) 

= 

1 .30003133E+0 1 

G(  9)  = 

9890=  G(  10)  = 
2.23737513E-07 

9.73213937E+03 

G( 

11) 

= 

-5.34933992E-09 

G<:  12)  = 

9100=  G(  13)  = 
9110=*EOR 

5.73437742E-09 

G< 

14) 

s 

-2 . 43 134724E-0S 

G( 

9120=1  CSA  NOS/BE  L534E 

L534 

CMR 1 

1/01/33 

9130=  15.38.50. LWMHUM4  FROM  /HU 

9140=  15. 33. 50. IP  00000320  WORDS  -  FILE  INPUT  ,  DC  04 
9150=  15. 38. 56. LWM.T35, 1 0100, CM  100000. T830229, MADDOX 
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Initial  Analysis  of  Results 

The  5  thruster  technologies  were  input  to  SUMT  one  by  one  and 
optimized  against  the  NavStar  GPS  mission  model .  Figures  9  and  lb  show 
a  compilation  of  the  results  o f  these  35  runs. 

The  output  o-f  the  runs  is  in  terms  of  a  minimized  mass.  Since 
basic  vehicle  mass  is  not  changing  as  greatly  as  power  supply  ana 
propellant  mass,  these  are  the  only  masses  minimized  by  the  objective 
function  in  initial  runs.  With  the  final  X  values,  the  transfer  times 
must  be  calculated.  A  programmable  calculator  handled  this  task  well, 
and  it  should  not  be  difficult  for  users  to  do  the  same  or  to  write  a 
FORTRAN  code  for  a  larger  machine.  Examples  of  how  this  transfer  time 
was  calculated  appear  in  Figure  11.  Follow-on  work  mav  be  able  to 
include  the  entire  methodology  in  a  single  code  which  calls  on  SUMT  and 
QGERT  as  subroutines,  thus  eliminating  the  need  for  separate 
calculations. 
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V/fcM&Les: 


Gea/eeal  Equat/okI: 

=  -transfer 

£"  =  Mp(K6)  __  Mp^e.  _  Mpiclsp  tfmt  (**c) 

rb  (*6/s)  ~  **T  ~  3*T 


(Mfl+Mp2)j<.T*P 

$*T 


Humier  fff 


Thrusitrs 


ft  hfcktr- 
rriiF 


Hn-»z£t 


rrv* 


fp  return. 


Example  •  908X6  p*pJo~fj  3-6*4  $/>** 


\L  s  eJ 

Vc.bc.tty 


(Mpi+Mp2)3cTSp  _  (*3  +  xs)  jU*) 

?*T  “  S  ft  T 


Z  - 
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Tsf  =  Sfcrjftt, 

I»f****- 

T  -  Thrust 
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J  crusfuuE. 


reJn+i-nrAfp^^. 


iT: 
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FIGURE  II.  PerERM/AJ/A/G,  TfAttSFEE  T/M& 


Note  that  the  transfer  time  calculations  were  made  for  point 
designs.  That  is,  when  SUMT  finished  optimizing  the  EGTV  for  a  given 
payload  mass,  the  transfer  time  was  calculated  based  on  tnose  final  X 
values.  This  is  irrespective  of  the  reduction  in  transfer  time 
possible  with  off-optimum  power  input  and  mass  change.  In  all  cases  it 
is  possible  to  reduce  these  transfer  times  by  increasing  the  thrust  or 
the  number  of  thrusters,  even  though  the  mass  increases  as  a  result  of 
using  more  power.  This  is  illustrated  in  the  sensitivity  analysis. 
Figure  12. 

Referring  back,  now,  to  Figures  ?  and  10,  it  can  oe  seen  that 
transfer  times  are  out  of  reason  for  the  Extended  Performance  thruster 
with  simplified  PPU  and  for  the  Argon  thruster.  The  reason  is  that 
minimum  mass  occurs  at  the  lower  limit  of  input  power.  Ms  input  power 
is  increased  above  this  lower  limit,  mass  of  the  power  supply  rises 
rapidly  .just  as  for  the  other  thrusters,  but  the  mass  of  the  propellant 
does  not  decrease,  but  increases  slightly  to  offset  the  increasing 
power  supply  mass.  Thus,  while  a  crossover  point  does  occur  between 
curves  in  the  feasible  region,  the  sum  of  the  power  supply  mass  and 
propellant  mass  increases  monotonically  throughout  the  feasible  region. 
Therefore,  the  SUMT  program  is  correct  in  finding  the  minimum  mass  at 
the  lower  limit  of  input  power,  2.0  KW,  but  the  result  is  not  very 
useful  when  considered  alone.  In  this  case,  the  decision-maker  must 
specify  a  desired  minimum  transfer  time  and  work  back  through  the  set 
of  equations  to  obtain  the  input  power  and  other  parameters  of  the 
EGTU.  Minimums  within  the  limits  of  the  feasible  region  do  occur  for 
the  other  three  thrusters  and  the  minimum  mass  and  resulting  X  vectors 
are  found  by  SUMT. 
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Basel  OH  JO -cm  3-G8/L  Thauste/l,  7000 7C6 
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xi-w 

3.700 

3.720 

3.  792. 

7.oo 

2o.  oo 

XS-Hn 

*733.33 

434.80 

452,20 

378.46 

301.31 

F -Mass 
OBJ  fed 

19 02.  7 

1902.72 

1902.56 

2 085.02 

3775.72 

Theust 

FEE 

EM6/M6. 

. 201237 

.202473 

.204772 

.206  2/0 

.  72o34o 

i y  ot*r 
w  (om) 

/44J7 

7  43, 77 

742.35 

J03.  38 

65.25 

fy.  fctur*. 

Cl  (ant) 

Mto 

74.68 

74.25' 

62,14 

49.47 

KouHb 

'c  r*» 

(Mis) 

212.97 

2/2.43 

2/6.  to 

765.5/ 

!  14. 72 

i 

hJoTE :  A8ev£  08  BELOW  730  0FT/A4mM j  F  (jCiA*m0*  MA*$  XdCAEAses. 

MAtti  TDEA  t  As  JM/ttr  TacAeasec.  Tsr  Aut>  thFuit  Jwcaease^ 

Aut>  TFAMSfEA  T/ME  STEAS/cy  ICctfASES.  LESA/rE  <SEEATF£. 

mass  of  7X0  A*#e<.  SifSAty. 


FIGURE  12.  TRANS  EE*.  71  HE  SF/ZS/T/y/ry 


Concentrating  now  on  the  final  minimized  mass  (F) ,  it  can  be  seen 
that  each  category  of  thruster  does,  in  fact,  require  a  different  mix 
of  prooellant  vs  power  supply  mass  to  deploy  the  same  set  of  GPS 
satellites.  As  discussed  in  the  previous  paragraph,  without  a  tradeoff 
of  power  supply  mass  and  propellant  mass,  both  Argon  and  Extended 
Performance  thrusters  have  minimum  mass  values  at  the  lower  limits  of 
mout  power.  Consequently,  they  have  the  smallest  objective  function 
mass  values  (F) ,  respectively.  The  next  lowest  is  Baseline,  followed 
by  Ring-Cusp  3-Grid  Xenon  and  3-Grid  rig.  With  its  high  thrust  to  power 
ratio,  the  Ring-Cusp  3-Grid  Xenon  has,  by  far,  the  snortest  trio  times. 
In  order  to  compare  the  5  thrusters  at  the  same  input  power,  tne  last 
table  in  Figure  10  shows  non-optimized  calculations  of  each  thruster 
operating  at  4.9  KW.  Each  is  carrying  a  2724  Kg  payload  representing  3 
GPS  satellites.  Each  has  the  same  power  supply  mass.  Each  deploys 
from  a  200  KM  Shuttle  orbit  to  GEO  and  returns  empty.  The  3  GPS 
satellites  are  assumed  to  be  deployed  in  the  same  orbit  plane  with  eacn 
satellite  using  its  ACS/RCS  to  position  itself  within  tne  orbit  plane. 
CFol low-on  studies  may  address  the  possibility  of  deploying  to  other 
orbit  planes  using  the  EOTV,  as  this  seems  quite  feasible  for  EP3.  hs 
the  figure  shows,  the  Ring-Cusp  Xenon  thruster  ranks  number  1  in 
shortest  transfer  time,  and  number  2  in  least  mass  to  accomplish  tne 
mission.  Its  relative  position  regarding  transfer  time  will  hold 
across  the  spectrum  of  payload  masses.  Its  relative  position  regarding 
minimized  mass  will  trade,  however,  with  at  least  one  other  thruster. 

The  question  arises,  "given  these  outputs  as  summarized  on  Figures 
9  and  10,  how  does  the  user  pick  an  optimal  EOTV  configuration!5"  If 
there  were  no  clear  cut  choice,  as  there  is  in  this  preliminary  set  of 
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runs ,  s  mini  mum  m-a  s  s  within  an  -acceptable  transfer  time  would  oe 


chosen.  For  instance,  say  3  GPS  satellites  form  the  payload  tor  which 


an  optimim  EOT'?  is  desired.  Optimized  final  masses  (F)  for  each 


category  thruster  are:  Baseline.  1635  Kg;  3-Grid,  2342  Kg;  Extd  Pert., 


1328;  Argon,  1983  Kg;  and  Ring-Cusp,  131'?.  These  represent  combined 


power  supply  and  propellant  masses  for  the  optimized  EOTV.  In  the  same 


order,  round-trip  transfer  times  are  370,  254,  484,  1880,  and  133  days, 


If  the  acceptable  transfer  time  cut-off  is  300  days,  then  2  choices 


remain;  3-grid  Hg  and  Ring-Cusp  Xe.  Between  these,  Ring-Cusp  uses 


less  mass.  Thus,  the  final  choice  for  an  optimized  EOTV  to  do  the 


mission  is  a  Ring-cusp  3-Grid  Xenon  Ion  Thruster  driven  at  4.312  KW 


with  a  solar  array  unit  of  748.293  Kg  and  using  341.53  Kg  of  propellant 


out  to  GEO  and  237.47  Kg  of  propellant  for  the  return  empty. 


In  this  case,  the  choice  was  obvious  because  Ring-Cusp  accomplishes 


the  same  mission  with  significantly  less  transfer  time  and  with  less 


propellant  and  power  supply  mass.  This  might  not  be  true  for  other 


choices  of  electric  thruster  technologies  which  a  user  might  wish  to 


evaluate.  Hence  the  need  for  examining  both  mass  and  transfer  time  in 


the  methodology. 


The  third  research  question  posed  in  Chapter  3  was,  "Does  one 


thruster  technology  clearly  outperform  all  others  for  each  mission?" 


It  appears  that  the  Ring-Cusp  design  is  a  clear  winner.  This  is 


consistent  with  the  fact  that  it  is  currently  being  investigated  as  a 


significant  improvement  to  the  baseline  38-cm  ion  thruster.  Given 


these  initial  results,  the  first  3  research  questions  have  been 


answered  and  the  first  three  objectives  of  the  thesis  have  now  been 
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CHAPTER  V .  Flv-Qff  Using  Q6ERT  Simulations 


Conceptual  Model 


The  -following  discussion  and  diagram,  Figure  13,  o-f  the  genera) 


conceptual  mode)  are  to  assist  in  understanding  the  -four  QGERT 
networks.  These  networks  model  the  -four  vehicles,  E0TU,  RBPU,  IUS,  and 
CENTAUR,  -for  the  Mf1y-off.“ 

Graphical  Evaluation  and  Review  Technique  (GERT)  witn  Queueing 
system  capability  (Q)  ,  or  QGERT,  is  a  modularized  and  easy  to 
■formulate  simulation  language  (37).  Follow-on  users  should  -find  QGERT 
easy  to  modi-fy  and  work  with  when  changing  the  models.  As  mentioned 
be-fore  in  Chapter  3,  the  strength  o-f  the  simulation  is  the  capability 


25 


i 


to  introduce  uncertainty.  Also,  the  capability  to  vary  inputs  to 


determine  long  term  (length  of  the  simulation)  effects  is  a  strength. 


In  order  to  show  the  rationale  for  each  model  and  the  general 


construction,  refer  to  Figure  14.  This  general  concept  applies  to  all 


four  networks.  The  Shuttle  turn-around  and  payload  integration  time 


are  external  to  the  system  modeled,  but  do  impact  the  Shuttle  launch 


rate.  The  model  begins  with  the  arrival  of  each  Shuttle  in  LEO 


carrying  one  or  more  payloads.  The  next  module,  payload  manifest,  is 


the  mission  model.  This  has  been  determined  previously  by  the  user. 


For  the  runs  made  by  this  author,  payloads  were  chosen 


probabi I istical 1 y  by  category,  realizing  that  not  all  payloads  needed 


transporting  to  higher  orbit.  The  next  module  represents  picking  an 


available  OTV  and  docking  prior  to  transfer.  This  module  does  not 


apply  to  expendable  upper  stages  (IUS  and  CENTAUR).  The  orbit  transfer 
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is  the  next  module  and  is  represented  primarily  by  the  time  delay 
required  -tor  the  operation.  Next  is  the  arrival  at  GEO  where  QTVs 
separate  and  return  via  the  QTO  return  module  to  a  queue  in  lEO. 
Statistics  are  collected  on  the  payloads  delivered  to  GEO  and  this 
mission  is  complete.  The  process  then  continues  with  the  arrival  of 
the  next  Shuttle  and  cargo.  With  this  general  conceptual  model  in 
mind,  understanding  the  individualized  networks  which  -Follow  should  be 
easier . 

Since  it  is  an  important  input  assumption,  the  mission  model 
applied  to  each  of  the  tour  vehicles  will  be  brie-fly  discussed.  For 
the  initial  runs,  a  mission  model  was  devised  which  used  the  masses  tor 
GPS  (Chapter  4)  but  accounted  tor  the  -fact  that  not  all  payloads  will 
be  candidates  -for  orbit  transfer  using  OTVs  or  IUS  or  CENTAUR.  Based 
on  numerous  mission  model  projections  -from  the  literature,  the 
following  assumptions  are  made: 

1.  65'/.  of  the  missions  brought  to  LEO  will  not  be  candidates  for 
using  the  fly-off  vehicles.  Reasons  include:  payload  is  to  remain  in 
LEO,  Spacelab  sortie,  PAM  is  being  used. 

2.  10X  of  the  missions  brought  to  LEO  aboard  Shuttle  carry  some 
other  LEO  payload  and  one  GPS  <983  KG)  or  a  similar  mass  satellite 
needing  transport  to  GEO. 

3.  20/i  of  the  missions  are  in  the  mass  category  of  3  GPS 
satellites  (2724  KG)  or  4088  lbs. (pushing  IUS  limits). 

4.  5X.  of  the  missions  are  heavy,  in  the  category  of  6  GPSs  (5443 
KG),  which  pushes  the  limit  on  an  updated  CENTAUR. 

Again,  the  user  is  free  to  use  another  appropriate  mission  model, 
but  this  one  had  the  nice  feature  of  “dual  representation."  That  is, 
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while  GPS  numbers  are  used  so  that  results  can  be  evaluated  tor  this 
system,  the  mass  categories  can  also  represent  other  satellites  just  as 
easily.  Examples  within  the  988  KG  range:  Landsat,  Nimbus,  2  or  3 
GOES,  Comstar,  RCA  Satcom,  Galaxy,  DSCS,  DMSP,  and  SOS.  Examples 
within  the  2724  Kg  range:  Newer  FI tSatCom,  TDRSS,  ERBS,  and 
combinations  or  multiple  satellites  in  the  908  Kg  category.  The  5443 
Kg  category  might  include  MILSTAR  and  multiples  o-f  the  988  or  2724  Kg 
groups.  Limitations  on  the  IUS  and  CENTAUR  are  one  payload  and  upper 
stage  per  launch,  currently.  Given  the  masses  involved  with  CENTAUR 
and  potential  payloads,  this  should  remain  a  good  assumption. 
Integration  that  must  take  place  on  the  ground  and  the  mass  o i  the  IUS 
also  make  it  likely  that  no  other  transportable  payloads  would  also  be 
aboard.  More  will  be  said  about  how  each  vehicle  handles  the  payload 
as  individual  networks  are  discussed. 


Electric  OTV  Network 

Figure  14  shows  the  QGERT  network  For  the  EQTV.  The  Shuttle 
arrival  rate  is  normally  (distributed  about  29.27  days  tor  13  launches 
per  year.  The  standard  deviation  of  this  arrival  rate  was  chosen  as  19 
days  to  introduce  more  variability.  This  will  likely  be  true  tor 
Shuttle  launches  well  into  the  future,  given  the  large  number  of 
variables  during  turn-around.  No  delay  exists  between  arrival  node  1 
and  payload  determination  node  4.  As  stated  previously,  35*/  of  the 
payloads  are  not  candidates  for  OTVs/ 1  US/'  CENTAUR  and  go  to  node  29, 
where  the  number  of  sorties  is  counted.  Then  19/,  29/,  and  5/  go  to 
nodes  5,  6  and  7,  respectively.  At  this  point,  the  payload  has  been 
selected.  Its  mass  is  assigned  in  the  appropriate  node.  This  mass  is 
later  used  to  determine  the  appropriate  transfer  time  in  the  QGERT  User 
Function  (UF)  .  Node  5  is  the  payload  scenario  where  one  GPS  or  similar 
mass  satellite  is  to  be  transported  to  GEO  and  the  remaining  cargo  bay 
space  is  used  for  LEO  satellites,  experiment  packages,  and  the  like. 
The  2724  Kg  node  in  the  EOTV  case  can  represent  actually  two  2724 
satellites  or  6  satellites  (this  is  not  possible  with  IUS  and  CENTAUR 
due  to  space  and  mass  limitations  mentioned  previously.)  Thus,  two 
activity  branches  join  nodes  6  and  19.  Node  19  is  the  queue  for 
satellites  ready  to  dock  with  EOTVs.  One  day  delay  is  built  into  the 
model  here  for  this  operation.  It  is  intended  that  Shuttle  and  crew  be 
present  for  monitoring  and  assistance  during  docking.  Node  11  is  the 
assembly  node  wich  selects  a  payload  and  an  EOTV  from  respective  queues 
and  begins  the  transfer.  At  this  point  the  QGERT  program  calls  UF  1 
and  assigns  the  appropriate  transfer  time  to  this  activity.  Transfer 
time  has  already  been  calculated  for  each  respective  payload  and  is 


input  in  number  of  days  to  the  UF  as  vanaole  Tl.  Arrival  at  GEG  is 
signified  by  node  12  and  apropriate  statistics  are  collected  by  nodes 
15,  16,  18,  19.  Node  13  is  the  return  path  tor  the  EOTV.  UF  2  is  the 
return  transfer  time,  also  previously  determined  and  input  to  UF  as  T2. 

Besides  the  uncertainty  modeled  in  the  Shuttle  arrival  rate  and  the 
payload  manifest,  uncertainty  is  also  included  as  a  reliability  figure 
for  the  EOTV.  The  figure  assumed,  .990,  arises  from  the  parallel 
redundancy  of  the  8  thrusters  and  associated  PPUs.  Reliability  for 
parallel  components  is  given  by: 

Rel <R)  =  1  -  <1  -  R)h  (4) 

Rel(R)  is  the  reliability  of  the  parallel  system,  R  is  the  reliability 
of  the  individual ,  identical  components,  and  n  is  tne  number  of 
identical  components.  Considering  .just  the  8-thruster,  4-BIM0D  subunit 
alone,  this  would  allow  a  reliability  as  low  as  .43  for  each 
thruster-PPU  combination,  if  the  subunit  were  to  be  .99  overall. 
However,  operation  on  only  one  thruster  prevents  total  mission  failure, 
since  the  vehicle  could  eventually  limp  back  to  LEO.  But  this  would  be 
far  from  desirable,  as  is  the  intent  of  equation  4.  Also,  the 
interface  module  and  avionics,  power  supply,  and  housekeeping 
functions/subunits  are  in  series  with  the  BIMODs,  so  that  this  thruster 
redundancy  is  mitigated. 

The  EOTV  thruster  technology  chosen  for  the  initial  runs  was  the 
optimum  picked  from  3UMT  resul ts:  the  Ring-Cusp  3-Grid  Ion  thruster 
operating  on  Xenon.  The  optimum  vehicle  configuration  was  for  the 
worst  case  payload,  5448  Kg.  Choosing  the  vehicle  optimized  for  tne 


heaviest  payload  means  that  lighter  payloads  will  be  delivered  raster 


than  if  vehicle  mass  had 

been  optimized 

for 

that 

1  l  gn  t  er 

pay  1 oad . 

Having  chosen  the  optimum 

vehicle  to  use 

for 

the 

f  1  y-of f  , 

transfer 

times  were  calculated  for 

each  payload 

and 

the 

return . 

For  the 

respective  payload  masses,  the  transfer  times  input  to  the  UF  Fortran 
IF  Statement  were  71.4,  114.4,  and  173.3  days.  The  return  time 
calculated  was  45.3  days.  This,  of  course,  does  not  vary  between 
payloads  since  each  vehicle  always  operates  at  a  constant  maximum 
thrust.  Solar  occultation  and  Van  Allen  degradation  have  not  been 
modeled  in  these  initial  runs.  Propellant  is  assumed  to  be  carried 
aboard  each  Shuttle  flight  which  bears  payloads  requiring  EGTV 
services.  Propellant  tanks  are  assumed  to  be  modularized  to  the  extent 
that  they  can  be  exchanged  via  the  Shuttle  manipulator  arm  or  with  EVA. 

Another  input  assumption  is  that  one  Shuttle  launch  is  required  for 
each  EQTV  to  deploy  the  optimum  11  vehicle  fleet.  More  will  be 
discussed  about  the  optimum  fleet  in  the  results  section  of  this 
chapter.  Input  card  listings  are  found  in  Appendix  IV  for  each  vehicle 
model  (Q6ERT  network) . 


Otv  Model  —  Qgert  network 


RBPO  Network 


Characteristics  of  the  RBPC  have  been  determined  primarily  from  a 
Boeing  study  (15)  and  -from  a  Systems  Engineering  study  (49),  with  some 
input  -from  other  references.  A  table  of  component  masses  and  more 
detail  as  to  this  chemical  OTU  will  be  found  in  Appendix  I.  Generally, 
it  is  to  be  a  reuseable  derivative  of  the  Centaur  using  upgraded  RL-10 
engines. 

Important  to  this  network  is  the  determination  that,  without 
aerooraking,  this  RBPV  must  have  2  Shuttle  missions  dedicated  to 
bringing  up  fuel  each  time  a  2724  or  5443  payload  is  to  be  transported 
to  GEO.  Even  the  908  Kg  mission  required  an  additional  Shuttle  flight 
with  fuel,  which  is  why  the  assumption  is  made  that  the  RBPC  would  not 
even  be  used  for  such  a  mission.  This  is  one  inflexibility  that  has 
been  accounted  for  when  comparing  against  the  other  3  vehicles. 

Referring  now  to  Figure  15,  the  network  is  similar  to  the  EOTV  in 
that  the  transfer  and  return  portions  are  essentially  the  same. 
However,  the  payload  module  or  mission  model  portion  of  the  network 
must  account  for  the  extra  Shuttle  missions  lost  to  refueling.  This 
has  been  modeled  by  nodes  21  and  22,  They  represent  Shuttle  launches 
carrying  a  specialized  refueling  pallet  with  pumps,  valving  and  tanks 
of  LOX  and  LH.  Mission  specialist  astronauts  who  are  qualified  for  the 
touchy  refueling  mission  must  also  be  aboard.  Both  nodes  29  and  5  oo 
not  utilize  the  RBPV,  so  connecting  activities  from  5  and  29  back  to 
node  1  only  complete  the  required  QGERT  arrival  scenario. 
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I  US  Network 


In  discussions  with  personnel  in  the  I  US  program,  it  seems  clear 
that  only  one  combination  of  payload  and  mated  IUS  may  fly  aboard  a 
given  Shuttle  mission.  Thus,  as  seen  in  Figure  16,  the  two  payloads 
which  can  be  handled  by  IUS,  90S  and  2724  Kg.,  are  taken  just  one  at  a 
time  by  a  single  mated  IUS.  Recall  that  the  EOTU  model  permitted  two 
2724  Kg  payloads  aboard  the  node  6  mission.  IUS  is  stretched  at 
present  to  transfer  2724  Kg.,  or  6000+  lbs.  But  it  is  assumed  that 
this  category  of  satellite  could  be  handled,  even  if  not  the  full  2724 
Kg.  The  5443  mission  is  beyond  the  capability  of  the  IUS  as  presently 
operated.  A  major  network  difference  from  the  OTU  may  be  noted.  There 
is  no  return  node  since  the  IUS  is  expendable.  The  IUS  reliability  is 
assumed  to  be  .965  and  node  14  serves  as  a  collection  node  for  those 
which  fail  in  the  29-year  simulation.  Characteristic  of  the  high 
thrust  chemical  Hohmann  transfer,  the  time  to  GEO  is  about  1/2  day  and 
is  not  significant  when  compared  to  the  large  transfer  times  with  EOTV. 
The  remainder  of  the  network  is  like  the  previous  ones  —  ie., 


primarily  for  statistics  collection 


FIGURE  16.  ZUS  AfOO£L  —  Qgert  network 


CENTAUR-G  Network 


The  CENTAUR-G  (CNTAR)  network  is  quite  similar  to  the  IuS  network 
in  arrival  rate,  lack  of  return  module,  and  a  transfer  time  of  .5  dav. 
(See  Figure  17.)  The  major  difference  is  that  the  CNTAR  would  not 
likely  be  mated  to  as  small  a  satellite  as  the  983  Kg.  mission.  Thus, 
it  is  assumed  CNTAR'S  large  size  is  a  point  o-f  in-flexibility  when 
rigidly  keeping  the  same  mission  set  for  each  of  the  flyoff  vehicles. 
It  is  also  stretching  the  current  published  CNTAR  capability  (11)  o-f 
10,080  lbs.  to  have  it  deliver  5443  Kg.  to  GEO.  But  it  is  assumed  for 
these  initial  runs  that  the  RL-18  engines  and  the  vehicle  will 
incorporate  design  upgrades  by  1999.  Reliability  is  modeled  as  .965. 


Fly-Off  Results  Summary 

Determining  the  optimal  fleet  size  for  the  EOTV  and  R6PV  was  the 
first  task  in  using  the  output  from  initial  runs.  Refer  to  Figure  13 
and  note  that  the  fleet  size  for  EOTV  was  varied  from  3  to  CO.  As  the 
numoer  of  parallel  servers  (EQTVs)  approached  11,  the  average  number  of 
satellites  waiting  in  LEO  (in  SATQ,  node  10)  for  transfer  dropped  to 
.0273.  The  average  waiting  time  in  the  queue  dropped  to  1.0839  davs. 
As  more  servers  (EOTVs)  are  added  above  11,  these  values  go  to  0  and  no 
satellites  ever  wait  for  transfer  —  an  EOTV  is  always  ready  in  the 
OTVQ,  node  9.  Below  11  servers,  the  wait  time  and  number  waiting 
begins  to  rise  exponentially.  Thus,  while  10  servers  might  be 
acceptable,  the  sensitivity  is  too  great  and  the  wait  time  climbs 
'apidlv  if  one  vehicle  fails.  Thus  the  optimum  number  of  EOTVs  for  the 


rhat  empties  the  satellite  oueue  was  11. 
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F I GU  RE  17.  Centaur- g  mooel  —  qgbrt  network 


Again,  referring  to  Figure  18,  it  is  noted  that  the  RBPV  curve  -tor 
average  number  of  days  wait  in  the  queue  tor  transfer  is  very  steep. 
Thus,  while  5  servers  (RBPMs)  would  be  acceptable,  the  sensitivity  is 
too  great.  That  is,  if  one  RBPU  fails,  the  queue  builds  very  rapidly. 
In  fact,  since  transfer  time  is  not  a  factor  for  this  chemical  system, 
the  reliability  figure  is  the  main  driver  in  setting  the  number  of 
RBPVs  for  the  fleet.  If  replacement  were  part  of  the  model,  this 
optimum  number  of  6  could  be  reduced.  Having  determined  the  optimum 
fleet  sizes  of  11  and  6  for  EOTV  and  RBPV,  respectively,  research 
question  # 6  has  been  answered  and  research  objective  #6  has  been  met. 

Table  3  summarizes  the  results  and  output  analysis  of  the  initial 
set  of  runs.  Approximately  40  runs  were  made  to  determine  optimum 
fleet  sizes.  But  after  that  determination  is  made,  only  one  run  for 
each  of  the  4  vehicles  is  needed.  Each  QGERT  run  simulates  20  years  of 
operating  the  fleet  and  also  repeats  the  20  year  simulation  50  times  to 
average  the  effects  of  uncertainty  and  randomness.  The  number  of 
repititions  required  had  been  determined  for  a  previous  study  using 
similar  models.  It  should  be  emphasized  that  the  results  were  based 
on  given  initial  input  assumptions.  These  assumptions  were  explained 
previously  in  the  vehicle  network  discussions.  It  should  also  be  made 
clear  that  the  results  are  somewhat  sensitive  to  the  input  assumptions, 
particularly  the  payload  mission  model  assumptions.  For  instance,  if 
more  satellites  were  assumed  to  require  use  of  transfer  vehicles,  the 
total  number  of  deliveries  would  go  up  and  more  EOTVs  would  be 
required,  though  perhaps  not  a  commensurate  increase  in  RBPUs.  If  a 
wider  variety  of  satellites  were  modeled  rather  than  using  the  "dual 
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TABLE  3.  Qgert  "fu/off"  results  summary 


representation"  or  representative  mass  assumption,  then,  again,  more 
EOTVs  would  be  required  and  perhaps  the  -fact  that  IUS,  CNTAR,  and  RBPV 
cannot  be  used  -for  the  whole  spectrum  of  masses  would  be  less 
noticable.  That  is,  the  gap  in  satellites  delivered  by  EOTV  and  the 
other  vehicles  would  be  narrowed. 

Table  3  indicates  that  each  vehicle  had  the  same  average  number  of 
Shuttle  launches  available  to  the  model.  A  percentage  of  those 
launches  carried  payloads  which  required  transport  to  GEO.  The  average 
requiring  transport  was  126.  One  exception  is  RBPV  which  nad  fewer 
payload-bearing  missions  due  to  the  requirement  for  dedicated  refueling 
missions.  Recall  from  the  RBPV  discussion  that  a  full  fuel  load  could 
not  be  carried  in  one  Shuttle  flight  for  either  the  2724  or  the  5448 
missions.  These  two  separate  Shuttle  launches  for  refueling  each  RBPV 
transfer  operation  significantly  reduces  the  number  of  satellites 
brought  to  LEO. 

Next,  note  the  number  of  satellites  delivered  to  GEO  bv  each 
vehicle  in  28  years.  This  is  an  average  number  of  satellites 
delivered,  since  it  varies  for  any  given  20-year  simulation  (based  on 
the  current  random  number  stream).  The  reason  fewer  are  delivered  than 
are  made  available  results  from  vehicle/payload  incompatibilities  and 
from  vehicle  failures.  It  can  be  seen,  however,  that  both  OTVs  launch 
more  satellites  than  seemingly  are  available.  This  is  because  node  o 
represents  one  2724  Kg  load  for  the  IUS  and  CENTAUR,  as  these  must  be 
launched  together  with  the  payload.  But  this  payload  bay  space  can 
used  for  another  satellite  or  for  3  GPSs  if  reuseable  OTVs  are 
utilized.  Thus,  node  6  represents  actually  two  times  2724  for  both 
OTVs.  Since  it  is  more  practical  to  launch  a  larger  mass  with  the 


RBPM,  both  2724  payloads  are  assumed  to  be  carried  on  one  vehicle, 
eliminating  the  need  +'or  refueling  two  RBPVs .  This  must  be  accounted 
for  by  the  user  in  the  analysis  since  the  model  does  not  handle  this 
necessary  quirk  in  the  RBPV  -formulation.  Again,  the  largest  difference 
in  average  number  of  satellites  delivered  stems  from  input  assumptions 
that  IUS,  CNTAR,  and  RBPV  have  limitations  as  to  the  payload  category 
which  each  can  handle,  either  by  design  limitation,  or  by  practicality 
limitation.  Example  of  the  latter  is  using  CNTAR  or  RBPU  for  small 
payloads  in  the  903  range.  Number  of  satellites  /  payloads  does  not 
give  a  complete  picture  of  capability,  since  some  payloads  are  much 
heavier.  Therefore,  the  next  entry  in  Table  3  is  the  total  mass  in  KG 
of  all  payloads  delivered  in  the  20-year  period,  averaged  for  50 
repititions  of  the  simulation.  Now,  CNTAR  outranks  IUS  as  would  be 
expected.  But  EGTV  still  has  the  lead. 

Launches  to  refuel  the  R8PV  definitely  reduce  the  capability  to 
deliver  as  many  satellites  to  orbit.  The  EOTV  requires  typically  1000 
Kg.  of  propellant  or  less.  This  is  carried  on  the  same  flight  as  the 
payload  in  a  modular  tank  assembly.  The  spent  tank  assembly  is 
returned  to  earth  for  filling. 

EQTV  can  be  seen  to  offer  potential  reductions  in  Shuttle  launches 
after  accounting  for  the  initial  11  launches  to  place  the  EOTV  fleet  in 
LEO.  Also,  it  potentially  can  deliver  more  satellites  to  GEO  in  a 
20-year  period  than  the  other  vehicles.  Research  questions  #4  and  #* 
have  been  answered,  and  objectives  #4  and  #5  have  been  met. 


CHAPTER  VI .  Overall  Result 


Assigning  Life-Cycle  Costs 

The  assignment  of  Life-Cycle  Costs  (LCCs)  to  each  of  the  vehicles 
being  compared  is  not  an  easy  nor  straight-forward  task.  It  is  made 
more  difficult  in  that  such  information  is  well  guarded  by  contractors. 
Often  it  is  just  as  sensitive  with  NASA  and  DOD.  Some  of  the  telephone 
conversations  on  the  subject  can  not  be  referenced.  But,  this  is 
perhaps  as  it  should  be  when  considering  the  legal  aspects  of 
contracting  and  when  considering  the  sensitivity  of  program  survival  to 
costs. 

Good  data  was  available  in  the  literature,  however,  and  forms  the 
primary  source  for  costing  the  models.  Some  studies  normalized  costs 
and  performed  sensitivity  analyses  to  determine  cost  effective 
directions  for  technology  development.  But,  the  intent  here  is  not  to 
again  seek  optimums,  but  to  attach  a  very  rough  estimate  of  LCCs  and 
compare  totals  for  each  vehicle. 

Table  4  contains  the  life-cycle  cost  summary  for  the  four  vehicles. 
IUS  and  CENTAUR  cost  figures  were  based  on  several  telephone 
conversations  (not  referenced  by  request)  and  on  several  literature 
sources.  Both  RBPV  and  EOTV  were  based  on  information  from  a 
combination  of  several  references.  The  EQTV  used  for  cost  estimating 
was  the  optimized  EOTV  modeled  in  the  QGERT  flyoff.  The  summary  table 
is  fairly  sel f-expi anatory .  Vehicle  costs  were  figured  and/or 
amortized  over  20  years  such  that  the  per-vehicle  cost  could  be 
multiplied  by  the  number  of  vehicles  required  during  the  flyoff.  The 
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number  of  Shuttle  launches  was  determined  by  QGERT  for  each  vehicle  , 
.vv  mission  set  combination.  Other  operational  costs  were  figured, 

particularly  for  the  reuseable  systems,  and  a  total,  29-year  LCC  was 
determined.  This  LCC  is  divided  by  the  total  Kgs.  of  payload  mass 
delivered  to  orbit  in  29  years  (from  QGERT)  to  obtain  a  ratio  of 
dollars  per  Kg  delivered  to  GEO.  This  cost  to  benefit  ratio  is  the 
basis  for  comparing  each  vehicle  over  the  life  cycle.  Based  on  this 
analysis  the  ranking  was:  EOTV,  CNThR,  REPJ ,  and  IU3. 


New  Missions  and  Enhancements 

New  missions  and  enhancements  have  been  suggested  by  the  results  of 
both  the  SUMT  and  QGERT  analyses.  Because  this  analysis  shows  a 
definite  cabability  to  do  orbit  transfer,  other  missions  involving  more 

r\ 

and  less  mass  were  examined  with  SUMT.  Results  aDpear  in  Tables  3-9. 

The  first  mission  investigated  was  the  Large  Space  Structure  vLSS; 
component  transfer  from  LEO  to  GEO.  Two  masses  were  chosen  approaching 
the  limit  of  one  Shuttle  load.  (Tables  5  and  a).  The  first  L3S 
payload  mass  of  20,009  Kg.  required  an  optimized  EOTV  with  higher 
specific  impulse  than  previous  payloads,  as  expected.  The  transfer 
time  out  to  GEO  is  nominal,  213  days,  and  the  return  is  quite  fast  for 
EP,  43  days,  since  24  thrusters  are  driving  a  light  load.  The  second 
payload  mass,  29,430  Kg.,  also  continued  the  trend,  requiring  yet 
higher  specific  impulse.  It  had  reasonable  transfer  times  of  223  days 
out  and  47.7  days  back  to  LEO  for  a  round  trip  of  approximately  275.7 
days. 

For  lighter  missions,  a  roving  intelligence  gathering  vehicle  with 
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TABLE  5 
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TABLE  6 
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TABLE  7 


Rover  Vehicle,  8  Rina-Cusp  XE  Thrusters. 
InterchanaeaDl e  Sensors.  580  K6 
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TABLE  7 


Rover  Vehicle,  8  Ring-Cusp  XE  Thrusters. 
Interchangeable  Sensors.  509  KG 
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TABLE  3 


Reoair/Re-furbish  Vehicle.  3  Rino-Cuso  XE  Thrusters. 
Interchangeable  Repair  Modules.  1808  KG 
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TABLE  9 


Free  Rover.  580  KG  Pavload.  8  Rino-Cuso  XE  Thrusters. 
Interchangeable  Sensors/Modules. 
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interchangeable  sensors  was  investigated.  It  is  seen  from  Table  7  that 
this  vehicle  would  be  able  to  travel  to  GEO  in  96  days  and  return  in  33 
days,  although  most  missions  would  probably  involve  closer  orbits.  The 
next  mission  investigated  (Table  3)  was  the  repair  /  visit  /  refurbish 
mission.  Again,  the  lighter  payload  allowed  more  propellant  to  be 
carried  and  the  vehicle  accomplished  transits  between  orbits  faster 
than  when  used  for  deploying  satellites.  The  last  mission  considered 
(See  Table  9)  involved  carrying  extra  fuel  for  several  LEO  -  GEO  - 
Return  trips.  Thus,  it  would  be  more  autonomous  and  posess 
multi-mission  capability.  The  first  leg,  carrying  the  most  propellant 
mass,  still  only  required  13S  days  transfer  time. 

It  was  not  necessary  to  use  QGERT  for  these  missions  since  the  IUS 
and  CENTAUR  were  not  competitors.  The  RBPV  has  such  a  large  fuel 
requirement  that  it  also  does  not  appear  to  be  a  contender  in  its 
present  parametric  form. 


It  appears  that  although  the  EOTV  can  offer  potentially  fewer  S/Kg 
for  delivering  payloads  to  orbit,  the  transfer  time  and  Van  Allen 
exposure  for  the  payload  owners  may  still  be  unacceptable  in  some 
cases.  An  all -EOTV  fleet  might  not  be  wise.  The  EOTV  is  a  strong 
contender  for  LSS  and  free  rover  type  missions,  since  chemical  vehicles 
use  enormous  amounts  of  fuel  in  the  former  case  and  have  a  greatly 
reduced  payload  fraction  in  the  latter  case.  For  those  payloads 
compatible  with  EOTV  transfer  times,  dollar  benefits  are  to  be  had. 
Mixed  fleets  were  not  specifically  addressed  by  the  methodology,  but 


could  be  modeled  by 

combining 

vehicles  and  missions 

in  one 

QGERT 

network.  From  the 

resul ts 

of 

both  SUMT  and  QGERT 

anal yses, 

and 

considering  the  low 

cost  of 

PAM- 

-D,  the  best  mix  of  upper 

stages  / 

QTYs 

/  rovers  appears  to 

be: 

1.  PAM-D  for  spinable  satellites. 

2.  CENTAUR-G  -for  heavier  rapid  transfer  payloads. 

3.  EOTV  -for  all  other  payloads,  using  several  for  rovers  — 
intelligence,  sensing,  refurbishment,  repair. 
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CHAPTER  VII 


.  Summary  /  Conclusions  /  Recommendations 


Summary 

The  need  which  gave  rise  to  this  thesis  is  the  need  to  enhance  or 
make  better  use  of  the  Space  Transportation  System  with  a  reuseable 
upper  stage  or  orbit  transfer  vehicle.  In  examining  the  mass  of  fuel 
required  to  operate  a  chemical  OTV,  it  quickly  becomes  apparent  that 
higher  specific  impulses  are  a  necessity.  Liquid  bi -propel  I  ant  engines 
have  pushed  the  theoretical  limits  of  specific  impulse  as  exemplified 
by  the  Space  Shuttle  Main  Engines.  But,  the  specific  impulse  needed 
for  practical  reuseability  in  near-earth  space  should  be  well  above 
that  which  is  possible  with  chemical  combustion  and  expansion. 

After  a  personal  visit  to  NASA-Lewis  Research  Center  in  summer, 
1983,  it  was  clear  that  a  great  many  refinements  have  taken  place  in 
electric  propulsion  technology  —  especially  electron -bombardment  ion 
thrusters.  These  thrusters  have  a  specific  impulse  normally  in  the 
2000  to  4000  sec  range  which  allows  mission  accomplishment  with  greatly 
reduced  propellant  mass  over  that  required  for  chemical  propulsion 
systems.  The  supporting  propulsion  module  with  avionics,  thermal 
control,  propellant  tanks  and  power  processing  has  been  developed  to  an 
advanced  state,  as  well  as  the  electric  thrusters.  This  system, 
developed  by  NASA,  is  envisioned  to  be  modular  with  two  to  ten  or  more 
thrusters  as  needed  to  cover  a  wide  range  of  thrust  requirements. 
Given  these  advantages,  it  seemed  a  good  candidate  for  Shuttle 
enhancements,  whether  as  an  upper  stage  or  as  a  repair/  refurbish/ 
retrieval  vehicle.  But,  several  issues  needed  to  oe  addressed:  Which 
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of  several  thruster  technologies  would  be  best  for  certain  missions? 
What  size  power  supply  would  be  required?  How  much  propellant  would  be 
required?  Numerous  studies  had  addressed  these  issues  in  one  fashion 
or  another,  but  none  had  performed  an  optimization  of  actual  thrusters 
/  prototype  technologies  followed  by  a  comparison  "fly-off“  against 
baseline  chemical  systems. 

The  decision-maker  investigating  Shuttle  enhancements  and  increased 
near-earth  capability  would  probably  like  to  find  the  optimal  electric 
system  for  a  required  mission  or  mission  set  and  "fly"  it  against  the 
current  upper  stages,  I US  and  CENTAUR-G,  for  comparison  of  performance 
and  cost.  Also  desirable  would  be  a  comparison  with  a  projected 
reuseable  chemical  system. 

Providing  the  decision-maker  with  this  type  of  information  has  been 
the  subject  of  this  thesis.  The  approach  has  been  to  first  develop  a 
method  for  parametrically  characterizing  existing  prototype  thrusters 
and  an  existing  prototype  solar  array  power  supply.  It  was  desired  to 
characterize  existing,  experimental  thrusters  rather  than  ideal , 
mathematical  projections  in  order  to  provide  the  decision-mader  with 
more  realistic,  conservative  data.  This  was  done  by  linearizing 
relationships  of  input  power  to  the  thruster  vs  ISP  and  vs  thrust  as 
obtained  from  measured  data  from  NASA-Lewis  Research  Center  and  Hughes 
Research  Laboratories.  A  very  key  point  in  this  thesis  is  that  these 
relationships  then  included  and  accounted  for  all  thruster  losses  and 
efficiencies!  The  rest  of  the  propulsion  subsystem  has  been  patterned 
closely  after  the  NASA  BIHOO  configuration.  The  primary  parameters 
then  input  to  an  optimization  program  (SUNT)  were:  power  supply 
specific  power;  thruster  input  power  vs  ISP  relationship  for  the  given 


A5J 


m 


thruster  technology;  and  vehicle  mass.  Relationships  from  the  rocket 
equation  were  also  incorporated  to  specify  propellant  mass  used  for  the 
mission.  The  mission  /  payload  was  characterized  by  two  parameters, 
Av,  or  velocity  change  required  for  the  orbit  transfer,  and  payload 


mass. 


The  specific  thrusters  chosen  for  the  optimizations  were  the 
baseline  NASA-Hughes  30-cm  J-Series  Hg  thruster,  the  30-cm  with  Argon 
propellant,  the  30-cm  with  3-grid  ion  optics,  the  extended  performance 
30-cm  with  simplified  PPU,  and  finally  the  Ring-Cusp  3-Grid  30-cm 
configuration  using  Xenon  propellant.  This  was  felt  to  offer  a  good 
spectrum  of  thruster  technologies.  (Foi low-on  users  may  select 
additional  ones).  The  specific  power  supply  was  the  NASA-Lockheed 
experimental  solar  array  with  a  total  system  specific  power  of  .052 
KW/KG.  The  mission  set  chosen  for  the  optimization  was  the  NavStar  GPS 
with  1-6  satellites  transported  at  a  time.  Carrying  six  satellites 
simultaneously  would  mean  deploying  1/3  of  the  GPS  constellation  to  an 
orbit  plane  and  using  either  the  OTV  or  the  satellite  RCS  to  achieve 
the  desired  position  in  the  orbit  plane. 

The  optimization  program  found  the  minimum  mass  to  accomplish  the 
given  set  of  missions  for  each  thruster  technology.  Minimizing  mass 
impacted  vehicle  cost,  launch  costs,  payload  capability,  and  transfer 
time.  In  the  program,  mission  performance  was  a  constraint  that  the 
vehicle  had  to  meet  —  it  had  to  provide  the  necessary  &v.  Since  the 
EQTVs  were  optimized  for  specific  payloads,  several  runs  had  to  be 
made.  But,  a  key  point  in  the  analysis  is  that  the  masses  in  the 
mission  model  could  correspond  not  only  to  GPS  satellites,  but  could 
represent  other  categories  of  payloads  as  well.  This  dual 
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representation  reduced  the  number  of  computer  runs  and  allowed  the 
analysis  to  indicate  EOTV  performance  over  a  wide  spectrum  of  potential 
payloads.  Once  outputs  from  the  35  runs  were  obtained,  transfer  times 
had  to  be  calculated.  Then,  with  both  optimized  mass  and  transfer  time 
as  criteria,  the  vehicle  offering  the  best  combination  was  picked.  Out 
of  the  thruster  technologies  analyzed  for  this  thesis,  the  clear  winner 
(most  optimum)  was  the  Ring-Cusp  3-Grid  Ion  thruster  operating  on  Xenon 
propellant. 

The  *fly-off*  simulated  operation  of  this  optimized  EOTV  for  a  20 
year  (user  selected)  period.  Also  simulated  were  the  IUS,  CENTAUR-6, 
and  a  Reuseable  Bi-Propellant  vehicle  based  on  CENTAUR  technologies 
with  the  RL-10  engine.  Long-term  (20-year)  operation  and  performance 
was  then  examined  without  having  to  build  and  launch  the  system.  The 
outputs  of  the  QGERT  ‘fly-off*  simulation  were  examined  and  analyzed 
for  these  results:  the  total  number  of  satellites  launched  over  the 
20-year  period;  the  potential  number  of  Shuttle  launches  saved;  the 
total  mass  (KG)  of  payload  placed  in  final  orbit;  and  the  number  of 
refueling  missions  required  for  the  RBPV. 

Using  these  results  from  the  fly-off,  the  assignment  of  rough 
Life-Cycle  Costs  (LCCs)  for  each  orbit  transfer  vehicle  system  was 
made.  These  costs  represented  R.D.T.6*  E.  not  yet  accomplished, 
acquisition  /  production  costs,  and  operation  costs.  Given  the  coarse 
assumptions  made  for  this  part  of  the  analysis,  the  results  of  LCC 
analysis  were  as  follows.  The  CENTAUR-G  upper  stage  had  the  lowest  LCC 
followed  by,  in  order:  EOTV;  RBPV;  and  IUS.  A  more  applicable 
figure  of  merit,  however,  was  the  ratio  of  dollars,  LCC,  per  KG  of 
payload  mass  delivered  to  final  orbit.  This  could  be  regarded  as  a 


-  90  - 


<r  *  -  .  «  a  *  ' 


cost  to  benefit  ratio.  For  this  figure  of  merit,  *LCC/KG  payload 
delivered,  the  ranking  was:  EOTV;  CENTAUR;  RBPV;  and  IUS.  CENTAUR 
and  EOTV  had  exchanged  first  place  ranking  because  EGTV  is  more 
flexible  in  delivering  a  wider  range  of  payload  mass. 

This  example  of  the  methodology  should  be  regarded  as  an  initial 
analysis  comparing  each  of  the  existing  and  proposed  vehicles.  The 
methodology  is  flexible  enough  that  different  mission  models  and 
accompanying  assumptions  may  be  incorporated  by  follow-on  users.  Other 
electric  propulsion  technologies  may  be  examined,  with  this  initial 
analysis  as  a  baseline.  Other  orbit  transfer  vehicles  may  be  compared. 
Thus,  besides  the  initial  results  obtained  showing  the  viability  of 
EQTVs,  the  methodology  and  algorithms  developed  should  prove  useful  to 
other  users,  planners,  and  decision-makers. 


The  results  have  been  presented  in  some  detail  in  the  Summary,  and 
it  should  be  helpful  to  the  reader  to  now  relate  these  results  to  the 
original  research  questions  and  objectives  which  were  delineated  in 
Chapter  3  of  the  thesis. 

The  first  research  question  posed  in  the  thesis  was,  "Which 
electric  thruster  technology  among  several  lab  prototypes  would 
optimize  an  OTV  in  terms  of  reduced  propellant  mass  and  reduced  power 
supply  mass  for  a  given  mission?"  Objective  #1,  choosing  thruster 
technologies  and  optimizing  them,  as  well  as  objectives  #2  and  #3,  were 
accomplished  using  SUMT  formulations.  Using  3UMT  results,  the  answer 


to  the  first  research  question  was  the  Ring-Cusp  3-0rid  Xenon  thruster. 
Round  trip  transfer  times  were  calculated  for  each  thruster  and  mission 
combination.  Questions  #2  and  443  dealt  with  transfer  times  and  with 


determining  a  single  best  thruster  technology  for  all  missions.  These 
were  answered  as  well.  Since  it  clearly  outperformed  all  others  in 
this  initial  study,  the  Ring-Cusp  thruster  was  consistently  the  optimum 
choice. 

The  next  three  questions  and  objectives  were  answered  and  fulfilled 
when  the  Q6ERT  models  were  developed  and  run  to  simulate  a  20-year 
“flyoff."  A  definite  reduction  in  Shuttle  launch  rate  was  possible 
when  using  resueable  EOTVs,  and  mitigated  for  the  chemical  OTV  due  to 
extra  fuel  launches.  The  optimal  fleet  size  determined  for  the  EOTV 
fleet  operating  without  IUS  or  CENTAUR  was  11.  The  optimal  chemical 
RBPV  fleet  size  was  6.  The  rough  estimates  of  LCCs  revealed  that  the 
vehicles  ranked  best  in  this  order  for  f/KG  delivered  to  GEO:  EOTV, 
CENTAUR-G,  RBPV,  and  IUS. 

The  last  research  question  and  objective,  dealing  with  new  mission 
possibilities,  both  were  accomplished  as  other  potential  missions  were 
examined.  The  SUMT  optimization  was  the  applicable  part  of  the 
methodology.  This  was  because  QGERT  comparison  runs  were  only  needed 
when  competition  existed  between  the  chemical  propulsion  vehicles  and 
the  EOTV.  EOTV  payload  ratios  were  clearly  superior  for  the  following 
cases.  Two  missions  delivering  LSS  components  to  GEO  were  input  for 
EOTV  optimization.  Also,  a  free-flying  rover  for  intelligence 
gathering,  remote  sensing,  and  satellite  repair  /  refurbishment  was 
optimized.  The  results  fell  within  the  feasible  region  and  show  that 
the  EOTV,  when  not  being  used  for  deployment  of  satellites,  could  be 
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used  -for  carrying  sensor  packages  over  “hot-spots",  replenisning 
modular  satellites  or  any  of  a  number  of  such  missions.  The  EGTV  has 
been  shown  to  be  able  to  carry  the  propellant  and  payload  •for  repeated 
trips  at  the  velocity  cnanges  required  by  these  missions. 

For  an  EQTV,  the  acceptability  of  the  long  transfer  times  ana 
radiation  exposure  in  the  Van  Allen  belts  must  still  be  assessed  by  the 
payload  owner.  Though  this  preliminary  study  indicates  cost 
effectiveness  for  an  all-EOTV  fleet,  a  mixed  fleet  of  CENTAURs  and 
EOTVs  might  be  more  effective,  given  that  perhaps  a  significant 
fraction  of  satellites  could  not  linger  in  the  Van-Alien  belts.  It 
should  be  noted,  though,  that  a  shielded  capsule  for  the  payload  might 
alleviate  some  of  the  radiation  and  would  be  feasible  given  the  results 
of  this  analysis.  It  was  shown  that  EOTVs  are  less  sensitive  to 
increases  in  payload  mass  than  are  chemical  propulsion  vehicles.  In  a 
mixed  fleet,  the  EGTVs  could  be  used  for  numerous  missions  when  not 
being  used  for  deployment.  As  mentioned  above,  orbiting  sensor 
packages  over  "hot-spots"  and  return,  repair  and  return,  and  retrieval 
are  capabilities  not  possible  with  expendable  stages.  Such  vehicles 
would  definitely  enhance  our  present  capabilities  in  near-earth  space 
as  well  as  create  new  capabilities  for  NASA,  private  industry,  and  DOD 
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Recommendations 

Given  the  results  o f  the  extensive  analyses  and  given  the  breadth 
of  the  literature  reviewed  -for  this  thesis,  three  sets  of 
recommendations  have  arisen.  The  first  set  has  come  from  seeing  the 
potential  uses  of  electric  primary  propulsion  for  expanding  U.S. 
capability  in  near-earth  space.  The  second  set  has  arisen  from  noting 
areas  of  the  methodology  that  can  be  improved  with  follow-on  work.  The 
third  recommends  how  actual  implementation  of  an  EQTV  capability  should 
begin . 

With  this  introduction,  the  following  three  sets  of  recommendations 
are  made  as  a  result  of  the  studies  performed  in  this  thesis: 

1.  A  phased  approach  should  be  adopted  to  bring  electric  primary 
propulsion  vehicles  into  general  usage. 

Phase  I  —  Develop  and  launch  an  on-orbit  prototype  3-thruster 
BIMOD  unit  with  the  following  specific  missions: 

(1)  Demonstrate  concept,  test  vehicle. 

<2)  Use  as  roving  sensor  platform  for  DOD.  Also  use  to  inspect 
malfunctioning  satellites. 

(3)  Measure  actual  Van-Alien  radiation  dosage  during  several 
trips  to  GEO  and  back. 

Phase  II  —  Place  two  more  vehicles  in  orbit  which  have 
appropriate  improvements  incorporated.  Primary  missions  should  involve 
intelligence  and  sensing. 

Phase  III  —  Place  a  small  fleet  of  vehicles  in  LEO  to 


complement  the  current  upper  stages  in  use.  Specific  missions  should 


(1)  Deployment  of  hardened  DOD  satellites,  including 
on-orbit  spares  above  GEO. 

(2)  Deployment  of  satellite  constellations  using  one  EGTV. 

(3)  Exchange  DOD  satellites  with  spares  on  an  irregular  basis 
to  extend  satellite  lifetime  and  to  thwart  unfriendly 
ASAT  planning  and  preplanning. 

(4)  Demonstrate  feasibility  of  disposing  of  nuclear  waste 
capsules  on  a  sun  intercept. 

(5)  Retrieve  satellites  for  inspection  or  refuroishment . 

(,6'j  Visit  NASA  multi-mission  modular  satellites  for  module 

replacement  or  exchange. 

(7)  Retrieve  spent  stages  for  possible  refurbishment  (CENTAUR) 
or  avionics  retrieval  (IUS).  Retrieve  dangerous  space 
debris. 

2.  These  follow-on  studies  are  recommended  and  may  be  regarded  a 
potential  thesis  topics: 

a.  A  more  complete  cost  study  using  the  QGERT  results. 

b.  Develop  an  optimization  scheme  for  Beginning-of-Life  (BOL) 
vs  End-of-Life  (EOL)  sizing  of  solar  arrays. 

c.  Take  the  entire  methodology  and  incorporate  it  into  a 
single  executive  computer  program,  thus  eliminating  hand 
calculations  and  the  many  separate  runs.  Develop  into 

a  management  information  system. 

d.  Develop  a  more  sophisticated  QGERT  model  with  more 
flexibility  for  handling  differing  mission  models  —  ie., 
more  modular  such  that  mission  changes  do  not  cause 


major  changes  to  the  model  network, 
e.  Continue  wording  with  SUMT  to  simultaneously  minimize  mass 
and  transfer  time  based  on  a  user's  weighting  of  the 
importance  of  each. 

3.  Given  the  potential  benefits  to  each,  both  NASA  and  00D  should 
jointly  fund  the  first  phase  suggested  above. 

This  third  suggestion  is  made  with  the  knowledge  that  any  new 
system  is  going  to  be  expensive,  especially  if  it  is  a  viable  space 
system.  However,  building  a  prototype  demonstrator  and  launching  it 
should  be  very  cost-effective.  If  program  funding  for  Phase  2  or  Phase 
3  is  slipped  or  delayed,  the  prototype  vehicle  could  still  be  used  as  a 
platform  for  testing  other  concepts  in  addition  to  its  Phase  I 
missions.  Additional  missions  might  include  demonstrating  modular 
repair  of  satellites,  retrieval  of  spent  satellites  or  space  debris, 
or,  in  fact,  any  of  the  missions  suggested  for  the  other  two  phases. 

These  recommendations  are  made  with  the  knowledge  that  electric 
propulsion  is  serving  very  well  at  this  writing  on  the  Navy  NWA 
satellite  program.  Pulsed  electric  thrusters  are  providing  secondary 
propulsion  for  stationkeeping  /  drag  make-up  for  this  highly  accurate 
navigation  satellite  system  (56) .  These  micropound  Pulsed  Plasma 
Thrusters,  though  not  envisioned  for  use  in  primary  orbit  transfer 
propulsion,  are  providing  greatly  improved  accuracy  for  the  NOVA  and 
the  Ballistic  Missle  Submarine  Fleet  it  serves.  An  improvement  in 
in-track  position  error  from  79  meters  in  a  0.3  day  period  to  less  than 
78  meters  in  a  6  day  period  has  been  achieved,  meaning  more  autonomy 


and  less  ephemeris  updating  (56) . 

The  results  of  this  thesis  show  that  besides  potential  economic 
advantages  to  employing  an  EOT')  -fleet,  the  increased  operational 
capabilities  suggested  by  this  system  would  be  the  greater  payoff. 
Instead  of  a  “push-the-button-and-watch-i t-go"  mode  of  operation,  the 
U.3.  could  move  toward  more  flexible  and  responsive  modes  of  operation 
in  near-earth  space.  Return  and  retrieve  features,  high  payload 
ratios,  and  large  velocity  change  increments  of  EOTVs  would  certainly 
enhance  the  present  capabilities  of  the  Space  Transportation  System. 


The  -following  drawings,  tables,  and  -figures  are  -from  the  references 
in  parentheses  on  each.  A  few  show  calculations  by  this  author. 
Figures  19  -  23  show  vehicle  configurations  for  all  but  IUS.  Tables  19 
and  11  show  masses  for  the  EOTV  propulsion  subsystem.  Figures  24  and 
25  and  Tables  12  -  16  show  the  thruster  data  provided  by  NASA-Lewis  and 
the  resulting  linearized  relationships.  The  ring-Cusp  data  was  already 
presented  in  Chapter  4. 

The  vehicle  configuration  for  the  Centaur-G,  as  shown  in  Figure  1? 
from  (11),  provides  an  idea  of  its  dimensions,  masses,  and  a  few 
subsystems.  The  IUS  is  not  included  because  it  is  operational  and  its 
configuration  generally  known.  Figure  20  shows  the  Boeing  Space-Based 
OTU  (15).  While  their  study  included  a  ballute  for  aerobraking,  note 
that  this  mass  is  deleted  for  the  oresent  configuration.  Upper 
atmospheric  heating  and  drag  are  not  yet  well  modeled  and  nearer  term 
technology  is  assumed.  Eventually,  aerobraking  technology  must  be 
developed  though,  if  manned  operations  are  to  be  realized  for  OTVs. 
The  masses  shown  were  used  as  a  rough  estimate  for  the  RBPV  and  as  a 
basis  for  fuel  requirements 

Figure  21  shows  t  , e  *1  layout  of  the  30-cm  Kaufman 
electron-bombardment  thrust  diagram  of  the  Ring-Cusp  thruster. 
The  iron  filing  map  shows  the  magnetic  field  line  enhancements.  Figure 
22  shows  the  BIMOD  unit  with  two  30-cm  thrusters.  This  was  assumed  to 
be  the  basic  thrust  subunit  for  the  EOTV.  The  vehicle  structure  above 
the  BIMOO  and  Interface  module  would  contain  solar  array  steering, 


Details  of 
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avionics,  housekeeping,  and  payload  interface  mechanisms, 
the  thruster  subunit  are  found  in  the  Design  Manual  (12).  Tables  il 
and  12  contain  the  mass  breakdowns  which  were  used  as  a  baseline  for 
eacn  of  the  five  thrusters  analyzed  with  SUIT.  Figure  23  shows  the 
modularity  of  the  BIMOD  engine  system  and  the  Interface  module.  Both 


the  EOTV  and  the  RBPV  have  been  represented  parametrically  for  the 


Centaur  G  weight  summary. 
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Centaur  Q 
Weight  (lb) 

Total  Contour  cargo  element  weight 

54,431 

Total  airborne  support  equipment 

7,462 

Spacecraft  airborne  support  equipment 

0 

Centaur  airborne  support  equipment 

7,462 

Total  vehicle  weight 

46,969 

Spacecraft  gross  weight* 

10,288 

Centaur  tanked  weight 

36.681 

Centaur  jettison  weight 

6,720 

Centaur  dry  weight 

6,163 

Centaur  residuals 

557 

Centair  expendables 

29,961 

PropeRants  (LH2  4  L02) 

29,707 

Engine  bums 

29,105 

Start  shutdown  4  vent 

602 

Hydrazine 

250 

Helium 
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Figure  23.  Thrust  Subsystem  Configuration  02) 
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APPENDIX  II. 


Use  of  Transfer  Curves 

Alfano  and  Wiesel  (1)  explicitly  solved  the  slow  timescale  optimal 
control  problem  f or  low  thrust,  minimum  time,  minimum  energy  orbit 
transfer.  Edelbaum  had  previously  solved  the  optimal  one-orbit  control 
problem  -for  the  first  time. 

Figure  27  represents  a  global  mapping  of  the  solution  space  for 
this  transfer.  The  mapping  is  in  semimajor  axis  —  inclination  space 
and  provides  explicit  total  velocity  change  requirements  for  any 
desired  transfer. 

Figure  26  provides  an  example  calculation  of  the  required  velocity 
change,  Av  for  208  Km  LEO  to  GEO  transfer  with  23.5*  inclination 
change.  Note  that  the  dynamics  are  independent  of  vehicle  specifics 
such  as  thrust,  payload  /  vehicle  mass,  and  specific  impulse. 
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Example  Calculat/om 

Qt'yiK  S  Any  /otoffr’USt  oJocfir/c.  frofuJs/ox.  OTV 

sr  L$7 TJjSb* 

LEO  —  2 O  a  Jem i 
GEO  =  &.L/0S2/  R§ 

=  2&>S°j  iHcJi**//* |  eJuotyo,  *r  .  7*f/2$  rrneJimms 

fte.uAe  27  y  y 

Find:  VtdoeJ+y  Tncrrmuf^  AV^  r*dfu/r*d  £*-  J4<l  -famt/ce- 

•Sblutitn:  p/nd/jy  cJtmmyt.  ix  semiMiJur-  axis  j  a  (H,) : 

^4**-«*  =  ^  +2#  bn  —  6S7S.  fife £*  =  d  PU%  ~  /  031367 

ahjku-  ■=■  C,L/062(  =  £.C/0£Zf  frcff 

d?Jo^  ~  S-0 

-  %  =  Mfiffbu* 

*1  F»/y  7te»;  iwt=^= 

-  y \tM/z.*Azs2  =  Sttoscjj**- 

Now,  Ftjttrt  27  i 

-  .7S  (x  At  not  ftmtjm  fitkH-  J>*rty  but  r4ti+r 

A__  _  ^  ^/>»t  "frmmtLrmt  -m'hVw  } 

X  =  “*  *  ' 


NotX,  Ccm.vtrb'hj  Jo  ^/fox.  i 

Z=-7s g?  » (ttnx'XX'sisa  3§) 

4V-( £w«3/»-— 

J3ctni/s  for-  fui/  usaje.  4«/ eyfi***Jt'on  of  Fj,  27,  f*i**d >'*  Ptf  0) 


FIGURE  26.  TfiAtfSAE*  CuAvE  CuculAT/pa/s 


APPENDIX  III.  SUMT  Example  Outputs 


Included  in  this  section  are  simply  Xerox  copies  of  some  of  the 
outputs  that  were  obtained  -from  the  SUMT  nonlinear  optimization  program 
runs  -for  the  EGTVs.  The  optimization  program  runs  consisted  of  the 
•five  thruster  technologies,  each  optimized  tor  the  GPS  mission  model. 
Thus,  each  EOTV  and  associated  thruster  technology  was  optimized  -for: 
a  983  kg  payload;  1316  kg;  2724  kg;  3632  kg;  4546  kg;  and  5443  kg 
payload.  Additionally,  a  1300  kg  payload  was  used  as  a  baseline.  Thus 
the  total  initial  number  of  runs  was  35.  To  avoid  bulk,  the  iterations 
which  SUMT  prints  out  between  initial  inputs  and  final  values  have  not 
been  included.  For  the  same  reason,  only  a  few  thrusters  and  Day loads 
are  represented  in  these  copies. 
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FIGURE  37  SuAt/H/ixy  Output  —  C Tutaup  Moe>£L. 


INPUT  CAPOS  »• 

SEN, 

«A000X,ELC0TV,S,2S, 1*33,13, 

SOU, 

1,0,1 ,0,N* 

NEC. 

2,i,i,p,H* 

ICS, 

3, 1 , 1 ,0,n» 

ST*, 

A,  1 ,  1  ,P,  B, 

0,4* 

ST*, 

5/182500,5 

1,0, A* 

9TA, 

A/LBSOOO ,  | 

1,0, A* 

ST*, 

7/LB |O«00, 

,1,0, A* 

9T«, 

8/9ATPAIL, 

, 1 ,0, A* 

SUE, 
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1,0. P, ,0,0.2, 11* 
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10/3ATQ, 0, 

,0,7,17,0,0,3,11 

SEL, 
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SN,CTC,H/1,,N, 10 

ST*, 
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1,0, A,, ,8/1* 

ST*, 

13/IETURN, 

, 1 ,P, *• 

ST*, 
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ST*, 
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1.0,1,40,3* 

ST*. 

I6/6E0-1, 1 

1,0,1, SO, IS*  . 

ST*, 
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,1,1,0, A* 

3TA, 
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, 1 ,0,8* 

*CT. 
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!,«,,, 2* 
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o.*s* 
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A,?,,,*,,* 

OS* 

ACT, 
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12.1*, ..23 
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1N/PINALST 
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S, 1 , CO, NOB 
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6.1, CO, 272 

• 
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7, 1 , CO, 5** 

• 

VAS, 
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•*•  EXECUTION  NRL  BE  ATTENPTEO  ••• 
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l,EN,M*Ul.l(J*,l<RPV,  1  1  ,  |<»,  19»3,  1S.0.  ,  7  3l'«.  50.S,o,2* 

SHU, 1,0, 1.0,** 

1 , 1  ,0,4* 

REG, 3,1,1 ,0,M» 

SI  A  ,  a ,  I  ,  1,P.M,9'I,«* 

S I A, 5/LR250U  ,  I  ,  |  ,  l>,  A* 

STA  ,  b/LHSODO  ,  t  ,  1  ,I>,A* 

ST A. 7/Lb  10000, I , 1 ,0, A*  . 

3TA,«/SAI7Alt,l,t,l),A* 

nut,q/uf ya,e,«,o,F,,o,o.?,  it* 

‘UIE,  l  0/S AT  0,0, 7. 0. F .  1  7,0,0.3,11* 

SEL,M/Ll**UP,ASM,CtC,K/1,.9,lO* 

ST  A  ,  1  ^/titO«A  ,  I,  l,t),A,,,H/l* 
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SI  A  ,  I  3/GEU-H,  I, 1,0,1), 40, 3* 

STA , 1 6/GtO” I , I, 1,0, 1, SO, IS* 
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ST  A  ,20/LEliPAiA,  1 ,  I  »D,  H» 

SIA,21/HEfUtLl,t,l,0,B«  - .  -  -  .  - . . -  -  - 

3tA.22/REU'EL2, 1, 1,0, B* 

AC  T , } , 4 , , , 2*  -  ■  -  — 

AC T , 4 , 20 , , , 30 , , 0 ,65* 

ACT, 4, 5,, ,7, ,0,10* 

ACT ,4,0, , , H, ,0,20* 

aL I , 4, / , , , 9, , 0  .  OS*  — — -  ...  . -  - - -  --  -  - - 

ACT, 20,  I, lit  1,1, SO* 

AC  T , S , i , NO ,1,51* 

AC!. 0,21, HU, 1,32* 
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ACT  .22, 1  .HU,  1  ,SS*  . . . ....  -  -  . — . 
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ACT, 1|, 12, CU, 0.5, tb/RUl VI XF* 
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ACI ,11 ,12,CU,0.5, 1B/PUTV3AE* 

ACT.ll, 12, Uf .  1.3I/0TV4XFEH*  .  .  .  . 

ACI  .  1 1 ,  12.IIE,  1 , 32/OTV5XEEH* 

ACT,1  1,1 2,UF ,1,33/UTVbXFEN* 
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FIN* 
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APPENDIX  V.  Other  Mission  Possibilities  /  Sensitivities 


The  first  two  -figures  contain  SUMT  results  tor  EOTVs  including  more 
realistic  calculations  of  power  available  to  the  thrusters. 
End-Of-Life  (EOL)  rather  than  Beginning-Of-Life  (BOL)  sizing  of  the 
solar  panels  has  been  included.  That  is,  the  effect  of  Yan-Ailen 
radiation  degradation  has  been  incorporated  such  that  the  full  power 
required  is  available  at  the  end  of  the  vehicle's  lifetime.  Provision 
has  also  been  made  for  more  avionics  and  housekeeping  power 
requirements  as  well  as  some  cabling  losses.  Further,  transfer  time 
calculations  include  Id’/,  occultation  during  earth  shadowing. 

The  next  five  figures  show  other  mission  possibilities  which  were 
discussed  in  Chapter  6,  but  with  the  earth  shadow  effects  included. 
These  missions  can  be  thought  of  as  representative  of  yet  other  rover, 
sensing,  repair  /  rescue  (satellite  only)  and  visit  type  missions  to 
enhance  operational  capabilities  in  near-earth  space. 
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TABLE  17. 

H390  KG  Pavload.  EQL  Sizino.  Qccu 1  tat  ion  Included 


9349=  ******************* 
****************** 

9350=  FINAL  VALUE  OF  F  =  4.22477445E+93 

9349= 

9370= 

9380=  FINAL  X  VALUES 
9390= 


************ 


9400=  X(  1)  = 
1. 10812144E+03 

3.03579551E+93 

XC 

2) 

= 

2 . 97974355E+02 

XC 

3) 

9410=  X<  4)  =  3. 1 1328282E+99 

9420= 

9430=  FINAL  CONSTRAINT  VALUES 
9440= 

X( 

5) 

7.32838 151E+82 

X  ( 

9450=  G(  1)  = 

2 . 37974855E+02 

2.53579551E+93 

GC 

2) 

= 

4.94420449E+03 

G< 

3) 

9440=  G(  4)  = 
8.89187834E+03 

9.70202315E+03 

GC 

5) 

= 

1 .098 12 144E+03 

GC 

6) 

9470=  GC  7)  = 
7.2233S151E+02 

2.4 1323282E+90 

GC 

3) 

= 

1. 43847 172E+91 

GC 

9) 

9430=  G(  10)  = 
2.59544322E-07 

9.24714135E+-93 

G( 

11) 

-1.72374753E-93 

GC 

12) 

9490=  G<  13)  = 

-2.15277394E-97 

GC 

14) 

= 

-4.73492527E-07 

GC 

9580=*EOR 


-135  - 


TABLE  18. 


2724  KG  Payload.  EQL  Sizing,  Occupation  Included 

9490=  H  H  !l  !(  H  H  )t  )t  H  H  )t  if  H  SI  H  H  H  H  H 

9590=  FINAL  VALUE  OF  F  =  4.57633125E+03 

9510= 

9520= 

9530=  FINAL  X  VALUES 
9540= 


9550=  X<  1)  = 
1.34059937E+03 

3.37346146E+83 

X( 

2) 

— 

3. 18957204E+02 

X< 

3) 

9530=  X<  4)  =  3.47330532E+00 

9570= 

9530=  FINAL  CONSTRAINT  VALUES 
9590= 

X< 

5) 

6.34 12375  lE+0  2 

X  ( 

9600=  G<  1)  = 

3 .03957204E+02 

2.37346146E+03 

GC 

2) 

= 

6.62153854E+03 

G  ( 

3) 

9610=  G(  4)  = 
3.659400  13E+03 

9.63104280E+03 

G( 

5) 

1 . 33059937E+03 

G( 

6) 

9628=  G<  7)  = 
6.74123751E+02 

2.97330562E+00 

G< 

3) 

= 

1.65266944E+01 

G( 

9) 

9630=  G<  10)  = 
-3 . 753 12084E-08 

9.31537625E+03 

G( 

11) 

= 

4.23460733E-09 

G< 

12) 

9640=  G<  13)  = 

-9.25061613E-07 

G< 

14) 

= 

-2 .66353 163E-06 

G( 

9650=*E0R 


TABLE  19 


LSS  2.  20.099  KG  Payload.  24  Ring-Cusp  XE  Thrusters. 
12  BIMQDs.  with  Qccultation 


10150=  HilHHHHHHHHHHHSiHHH 

10160=  FINAL  VALUE  OF  F  =  6.24333310E+03 

10170= 

10130= 

10190=  FINAL  X  VALUES 
10200= 


10210=  X(  1)  = 
2.73374196E+03 

6.06485577E+03 

X< 

2) 

— 

1 .2 1072537E+02 

X< 

3) 

— 

10228=  X<  4)  = 
10239= 

6.29577216E+00 

X( 

5) 

6 . 98900253E+02 

X( 

10240=  FINAL  CONSTRAINT  VALUES 

10250= 

10269=  G<  1)  = 

5 . 56435577E+03 

G( 

2) 

_ 

3.935 14423E+03 

G( 

3) 

1. 11072537E+82 
10270=  G<  4)  = 

9.37392746E+03 

G( 

5) 

— 

2.72374 196E+03 

G( 

6) 

s 

7.26625304E+03 
10230=  G<  7)  = 
5.98900253E+02 

5.79577216E+08 

G( 

3) 

= 

1 . 37042273E+0  1 

G( 

9) 

= 

10290=  G<  10)  = 
-4 . 57452325E-06 

9.39109975E+03 

G< 

11) 

= 

1.33833609E-07 

G( 

12) 

10300=  G(  13)  = 

-8 . 939350 66E-0  7 

G( 

14) 

=s 

-3.78521576E-06 

G( 

10318=*EOR 


TO  GEO  =  241.23  DAYS 
RETURN  =  53.50  DAYS 


ROUND 

TRIP  =  293.73  DAYS 


TABLE  20 


LS3  3.  29,480  KG  Payload,  32  Rina-Cuso  XE  Thrusters 


18  Bilious,  with  Occultation 


18300=  ****************  **************** 
****************** 

10310=  FINAL  VALUE  OF  F  =  8.81802384E+03 

10320= 

18330= 

10340=  FINAL  X  VALUES 
10350= 

10380=  X<  1)  =  8 . 23038055E+03  X<  2)  =  1 .2543729 1E+02  X<  3)  = 

3.30555340E+03 

10370=  XC  4)  =  8.5227391 1E+00  X<  5)  =  7.98432131E+02  X< 

10380= 

10390=  FINAL  CONSTRAINT  VALUES 
10400= 


1. 1543729 1E+02 
10420=  G<  4) 
8.19444880E+03 
10430=  G<  7) 
7 .38482 13 1E+02 


1.7 10984 17E-07 
18450=  G<  13) 
10480=*EOR 


5.73033055E+03 

G< 

2) 

= 

3.7191 1945E+03 

Gt 

3)  = 

9. 3745827 1E+ 03 

G( 

5) 

= 

3.79555340E+03 

G( 

8)  = 

8.0227391 1E+00 

G( 

3) 

= 

1 .34772809E+0  1 

G( 

9)  = 

9.2035 1737E+03 

G( 

11) 

- 

1. 833022 17E-09 

G( 

12)  = 

-3. 55357770 E-08 

G( 

14) 

s 

-2.93857849E-07 

G( 

TO  GEO  =  250.31  DAYS 
RETURN  =  52.49  DAYS 


ROUND 

TRIP  =  303.30  DAYS 


1 

\ 


TABLE  21 


Rover  Vehicle.  8  Rina-Cuso  XE  Thrusters 


Interchanaeabl e  Sensors.  590  KG.  Shadowin 


9220=  H  if  Hii  HHHHHHHX  HHUini  XH 

XXXXXXXXXXXXXXXXXX 

9230=  FINAL  VALUE  OF  F  =  1 .37461359E+03 

9240= 

9250= 

9260=  FINAL  X  VALUES 
9270= 


4.0 1332023E+02 


9300= 

9310=  FINAL  CCt 
9320= 

9330=  G(  1)  = 
6.84975642E+01 
9340=  G<  4)  = 
9.59886798E+03 
9350=  G(  7)  = 
3.3530 1052E+02 
9360=  G(  10)  = 
-4.79021764E-08 
9370=  G(  13)  = 
9330=XEOR 


3.95763998E+03 

X< 

2) 

= 

7.34975642E+0  1 

X< 

3) 

4.08187334E+00 

X  ( 

5) 

s 

3.4530 1052E+02 

X< 

5TRAINT  VALUES 

3.457S3990E+03 

G( 

2) 

= 

6.042310 10E+93 

G< 

3) 

9.92150244E+83 

G( 

5) 

= 

3.91332023E+92 

Q( 

6) 

3.58137334E+00 

G( 

8) 

= 

1 .59131267E+0  1 

G< 

9) 

9.65469895E+03 

G( 

11) 

= 

-2.61934474E-10 

G( 

12) 

-6.05959940E-98 

G< 

14) 

s 

-8.06412572E-93 

G( 

TO  GEO 

RETURN 

ROUND 

TRIP 


105.67  DAYS 
90.92  DAYS 

196.59  DAYS 


TABLE  22 


Reoair/Re-furbish  Vehicle.  3  Rino-Cuso  XE  Thrusters. 
Interchangeable  Repair  Modules.  1990  KG.  Shawdowino 


9550=  HHHHHHHHHHHHH  HH  H 
HHH3 ************ 

9540=  FINAL  VALUE  OF  F  =  1 .53783934E+03 

9570= 


9538= 

9590=  FINAL  X  VALUES 

9400= 

9410=  X(  1)  = 

4.44852544E+03 

4.41 148314E+02 
9420=  X<  4)  = 

4. 418584 12E+00 

9430= 

9440=  FINAL  CONSTRAINT  VALUES 

9450= 

9440=  G<  1)  = 

3.94852544E+03 

7.38139255E+0  1 
9470=  G<  4)  = 

9.911 13 107E+03 

9.55835149E+03 
9480=  G(  7)  = 

4. 1 18584 12E+00 

3.74139445E+02 
9490=  G(  10)  = 

9.41384035E+03 

4.4444871  IE-08 
9700=  G<  13)  = 

-1. 425035 13E-03 

97 10=XEOR 

X< 

2) 

— 

8.38139255E+0 1 

X( 

3) 

X< 

5) 

S 

3.34 139445E+02 

X( 

G( 

2) 

= 

5.53147434E+03 

G( 

3) 

G( 

5) 

= 

4.31 148314E+02 

G( 

4) 

G< 

8) 

= 

1.533 14 159E+01 

6< 

9) 

G< 

11) 

= 

4.33447074E-09 

G( 

12) 

G< 

14) 

, 

-2. 474553 18E-08 

G< 

TO  GEO  =  114.22  DAYS 
RETURN  *  181.73  DAYS 
ROUND 

TRIP  =  217.95  DAYS 
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TABLE  23 


Free  Rover.  598  KG  Payload.  8  Rino-Cuso  XE  Thruster 


Interchanaeabl e  Sensors/Modules 


1680  KG  Extra  Proeellant.  Shadowin 


9138=  H  H  H  H  H  H  H  H  H  H  H  H  i(  H  H  H 

9140=  FINAL  VALUE  OF  F  =  1 .61 174152E+03 

9158= 

9160= 

9170=  FINAL  X  VALUES 
9180= 

9190=  X<  1)  =  4 . 65776684E+83  X<  2)  =  9.26425162E+8 1  X(  3) 

5.65150266E+02 

9200=  X(  4)  =  4.31741084E+00  X(  5)  =  3.05451 120E+02  X( 

9210= 

9220=  FINAL  CONSTRAINT  VALUES 
9230= 

9240=  G<  1)  =  4. 15776684E+03  G<  2)  =  5. 342233 16E+03  G(  3) 

8.26425 162E+ 01 

9250=  G<  4)  =  9.98735748E+03  G<  5)  =  5.55150266E+02  G(  6) 

9 . 43484973E+03 

9260=  G<  7)  =  4.3174 1084E+00  G<  3)  =  1 .51825392E+8 1  G<  9) 

2.95451 120E+02 

9270=  G<  10)  =  9.69454888E+03  G(  11)  =  2.83726313E-03  G(  12) 

3.29460458E-09 

9280=  G(  13)  =  -2.33522240E-97  G<  14)  =  -6.70475S56E-87  G( 

9290=*EOR 


TO  GEO 
RETURN 


148.92  DAYS 
80.39  DAYS 


ROUND 

TRIP 


229.31  DAYS 


|Kn 
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